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Abstract — For LDPC codes operating over additive white 
Gaussian noise channels and decoded using message-passing 
decoders with limited precision, absorbing sets have been shown 
to be a key factor in error floor behavior. Focusing on this 
scenario, this paper introduces the cycle consistency matrix 
(CCM) as a powerful analytical tool for characterizing and 
avoiding absorbing sets in separable circulant-based (SCB) LDPC 
codes. SCB codes include a wide variety of regular LDPC codes 
such as array-based LDPC codes as well as many common 
quasi-cyclic codes. As a consequence of its cycle structure, each 
potential absorbing set in an SCB LDPC code has a CCM, and 
an absorbing set can be present in an SCB LDPC code only if 
the associated CCM has a nontrivial null space. 

CCM-based analysis can determine the multiplicity of an ab- 
sorbing set in an SCB code and CCM-based constructions avoid 
certain small absorbing sets completely. While these techniques 
can be applied to an SCB code of any rate, lower-rate SCB 
codes can usually avoid small absorbing sets because of their 
higher variable node degree. This paper focuses attention on the 
high-rate scenario in which the CCM constructions provide the 
most benefit. Simulation results demonstrate that under limited- 
precision decoding the new codes have steeper error-floor slopes 
and can provide one order of magnitude of improvement in the 
low FER region. 

Index Terms — Cyclic consistency matrix, absorbing set, error 
floor, quasi-cyclic LDPC codes. 

I. Introduction 

Low-density parity-check (LDPC) codes were introduced 
by Gallager and are well-known for approaching capacity with 
iterative decoding. However, in the low FER region a flattening 
of the frame error rate (FER) curve called the error floor 
usually occurs for LDPC codes with moderate block lengths 
and high code rates. This performance degradation is due at 
least in part to the sub-optimality of message passing in a 
graph with cycles. 

This paper focuses on one class of regular LDPC codes: 
separable, circulant-based (SCB) codes. SCB codes include 
array-based LDPC codes as well as many common quasi- 
cyclic codes. SCB codes retain standard properties of quasi- 
cyclic LDPC codes such as girth, code structure, and com- 
patibility with existing high-throughput hardware implemen- 
tations. These codes might be applicable for high-throughput 
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data storage applications when message passing decoding is 
used with limited precision. 

This paper introduces the cycle consistency matrix (CCM) 
as a powerful analytical tool for characterizing and avoiding 
any graphical structure in SCB LDPC codes that contains a 
well-defined cycle structure. For LDPC codes operating over 
additive white Gaussian noise channels and decoded using 
message-passing decoders with limited precision, absorbing 
sets have been shown in (TJ, 12, |3| to be a key factor in error 
floor behavior. As a consequence of its cycle structure, each 
potential absorbing set in an SCB LDPC code has a CCM, 
and an absorbing set can be present in an SCB LDPC code 
only if the associated CCM is not full column-rank. 

For a specified circulant matrix, all SCB codes share a 
common mother matrix. For a specified absorbing set, the 
CCM approach yields with relative ease (because testing for 
a nontrivial nullspace is straightforward) a complete charac- 
terization (as to which codes have the absorbing set) for the 
entire family of SCB codes of a specified rate obtained by 
selecting a specified number of rows from a specified SCB 
mother matrix. The CCM-based analysis shows that quasi- 
cyclic code families described in 0), |5| and J6) include 
codes with good absorbing set spectra with a proper choice of 
parameters. This paper also identifies absorbing-set-spectrum 
equivalence classes within a family of SCB codes, significantly 
simplifying this characterization. 

Larger absorbing sets are more difficult to avoid, and our 
characterizations reveal that absorbing sets of large enough 
size occur in every code in the family. A specific code avoiding 
as many small absorbing sets as possible can be selected for 
further improvement. This improvement can be accomplished 
by existing construction techniques such as JH] described 
in Section [TTJ Alternatively, a CCM-based code shortening 
technique that we introduce in this paper can remove variable 
nodes so as to force the null space of the CCM to have zero 
dimension thereby avoiding the target absorbing set while still 
precluding the smaller absorbing sets already absent before 
shortening. 

Section [TT] reviews trapping sets, absorbing sets, and prior 
work on designing LDPC codes that avoid these structures. 



Section III introduces separable circulant-based (SCB) codes 



and the cycle consistency mahix (CCM). 

For SCB codes with column weights r = 5, Section IV 
considers selected-row (SR) SCB codes, which are constructed 
by selecting r = 5 appropriate row groups from the SCB 
mother matrix. This section identifies the CCMs for the 
smallest absorbing sets for r — 5 SR-SCB codes and uses 
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these CCMs to select rows that avoid all (4,8), (5,9), and 
(6,8) absorbing sets. Recognizing equivalence classes of SR- 
SCB codes that have the same absorbing sets greatly simplifies 
analysis and permits a complete characterization of r — 5 
SR-SCB codes with respect to the smallest absorbing sets. 



Section IV concludes with an example characterization of all 
r = 5 SR-SCB codes with respect to (4,8), (5,9), and (6,8) 
absorbing sets for code length n = 67 2 = 4489. 

Section [V] identifies the CCMs for the smallest absorbing 
sets for r = 4 SCB codes. These configurations are (6,4) 
absorbing sets. In this section, the ability to characterize an 
entire SR-SCB family is utilized to show that one of the 
three configurations is necessarily present in every r — 4 SR- 
SCB code regardless of p. The other two configurations can 
be avoided by proper row selection. This section concludes 
by showing how the CCM approach can guide the removal 
of specific columns from an SR-SCB parity-check matrix to 
eliminate the troublesome (6,4) configuration, producing a 
shortened SR-SCB code with no (6,4) absorbing sets. 

Section [Vljbriefly remarks on the r = 3 case, completing an 
analysis of small absorbing sets in SCB codes with 3 < r < 5. 
It is important to recognize that the CCM approach itself easily 
generalizes to any r. This paper focuses on 3 < r < 5 to 
provide illustrative examples and because r — 4 and r = 5 
are important for high-rate codes that are of interest to the 
authors. Section [VTIJprovides simulation results demonstrating 
performance improvement in the error floor position and slope 



obtained by the newly identified codes. Section VIII delivers 
the conclusions. 

Conference presentations of parts of this paper appear in 
HH, ED, G2, and ma. 

II. Trapping sets and absorbing sets 

Prior work indicates that certain sub-graphs called trapping 
sets dJ4), and, in particular, a subset of trapping sets called 
absorbing sets lfl5ll are a primary cause of the error floor 
in practical implementations. Absorbing sets are trapping sets 
that are stable under bit-flipping decoding. For LDPC codes 
operating over additive white Gaussian noise channels and 
decoded using message-passing decoders, some absorbing sets 
can be successfully overcome with sufficient precision. How- 
ever, in the practically important scenario of limited precision, 
absorbing sets have been shown in (TJ, 12), O to be a key 
factor in error floor behavior. 

We refer to the smallest absorbing sets (in terms of the 
number of variable nodes in the set) as dominant because 
an absorbing set is activated when all of its variable nodes 
are in error and this activation is more likely for a smaller 
number of variable nodes than a larger number of variable 
nodes. Our simulation results also show that when precision 
is sufficiently limited, error performance is dominated by the 
smallest absorbing sets in the cases that we studied. However, 
there may be situations in which a relatively large absorbing 
set structure with high multiplicity could be more important 
to performance than an extremely rare absorbing set that is, 
nonetheless, smaller. 

Since the performance degradation due to trapping sets 
(including absorbing sets) is partially due to the sub-optimality 



of the iterative decoding algorithm (especially under limited 
precision), one possible direction to mitigate the effect of 
trapping sets is to improve the decoding algorithm. More ef- 
fective message-passing algorithms ifTBI . ifTTl . better iteration 
averaging schemes [18], more efficient quantizations |[T9l . (3l 
and post-processing for absorbing sets [2| all can improve the 
error floor. 

A complementary direction to improving the error floor, 
and the focus of this paper, is to design a parity check 
matrix to avoid the trapping sets and absorbing sets. Numerous 
previous papers have taken this approach. The Approximate 
Cycle EMLjJ (ACE) algorithm |2U) or the Progressive Edge 
Growth (PEG) algorithm 11211 can be used for a column- 
by-column constructions. Algebraic methods for constructing 
LDPC codes also demonstrate good error floors as compared to 
randomly constructed codes. Results in 11221 . Il23ll . [24 1, |25|, 
and ll26l present some notable work in this area. 

Recent papers have proposed methods to improve the ab- 
sorbing set spectrum. Introducing additional check nodes E71l 
or increasing the girth ll28l eliminates small trapping sets 
for some codes. The algorithm in l26l and ll29l constructs 
quasi-cyclic codes from Latin squares so that the Tanner 
graph does not contain certain trapping sets. In ]7) §0, a 
deterministic construction of structured LDPC codes free of 
dominant trapping sets was developed by carefully swapping 
edges in the lifted graph of the original code to provably 
eliminate the presence of detrimental trapping sets. 

Several of these approaches have taken advantage of the 
cycle structure of the target objects. For example, the ACE 
algorithm avoids stopping sets (and incidentally trapping sets 
and absorbing sets) by ensuring that all small cycles have 
sufficient extrinsic message degree. Also, the cycle structure of 
trapping sets is utilized in Q JH] to inform the edge swapping 
in the lifted graph. In fl5l . Dolecek et al. used the cycle 
structure to study absorbing sets in array-based LDPC codes. 

Our work builds on fUl by capturing the relevant cycle 
structure of a specific absorbing set with a cycle consistency 
matrix and ultimately by whether the nullspace of that matrix 
is zero or not. A unique feature of the CCM approach is that 
it allows the characterization of an entire family of SCB codes 
rather than focus on constructing a single code by a column- 
by-column construction as in ACE or PEG or by a specific 
lifting as in Q QQ. ° n me otner hand, in this paper the CCM 
approach restricts attention to SCB codes while ACE, PEG, 
and the approach of Ql |8] can be applied in principle to 
any LDPC code. While the rest of the paper applies the CCM 
approach to absorbing sets, it can as well be applied to any 
trapping set or more generally any graphical object that is 
comprised of specific cycles connected in a specific way. 

III. SCB Codes and the Cycle Consistency Matrix 



Section III-A introduces separable, circulant-based (SCB) 
codes. Section III-B| introduces absorbing sets and then puts 
forth the main concept of this paper: the cycle consistency 



matrix (CCM). Section III-C uses graph theory to identify 



efficient CCMs that have the smallest possible number of rows. 

'EMD is the acronym for Extrinsic Message Degree. 
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S ection [IH-D | shows how the VN-graph approach simplifies the 
identification of efficient CCMs when the check node degree 
of the graphical structure is constrained to be two or less. 



Section III-E concludes the section with a theorem providing 
necessary and sufficient conditions for the existence of a given 
absorbing set in SCB codes in terms of its CCM. 

A. Separable, circulant-based LDPC codes 

Circulant-based LDPC codes are composed of circulant 
matrices and form a subset of (r, c) regular LDPC codes, 
where r is the variable-node degree and c is the check-node 
degree. Each circulant matrix is a p x p matrix, where p is 
a prime number. Parameters r and c are positive integers and 
are at most p. The structure of these codes is compatible with 
high-throughput hardware implementations (3j. 

The parity-check matrix of circulant-based LDPC codes can 
be described as follows: 



H 



p,f 



T f(o,o) 
r /(i,o) 

T /(2,0) 



r /(0,l) 



T /(2,l) 



r /(0,2) 
r /(1.2) 
T /(2.2) 



r /(0,c-l) 
r /(l,c-l) 
,/(2,c-l) 



where a is the following p x p circulant matrix: 
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1 



and f(i,j) is any function mapping the (row-index, column- 
index) pairs to the integers {0, . . . ,p — 1}. 

A column (or row) group is a column (or row) of cir- 
culant matrices. Each variable node has a label (j, k) with 
j € {0, c— 1} being the index of the corresponding column 
group and k E {0, ...,p — 1} identifying the specific column 
within the group. Similarly, each check node has a label 
where i G {0, r — 1} and I G {0, ...,p — 1}. 

This paper focuses on separable, circulant-based (SCB) 
codes, which are defined as follows: 

Definition 1 (Separable, Circulant-Based (SCB) Code): 
An SCB code is a circulant-based LDPC code with a 
parity-check matrix H^'f in which f(i,j) is separable, i.e., 
f(i,j) = 9r{i) ■ 9c(j) mod p. ■ 

Parity check matrices of SCB codes with a specified circu- 
lant matrix (i.e., a specified p) can be viewed as originating 
from a common SCB mother matrix with f m (i,j) = i-j 

mod p. The functions g r (i) and g c (j) effectively specify 
which rows and columns of the mother matrix are selected 
for the resultant SCB matrix. The ranges of g r (i) and g c (j) 
are both {0, . . . ,p — 1}. 

SCB codes include, for example, the constructions in |4j, 
0, and |6).The girth of all SCB codes is guaranteed to be at 
least 6 by the SCB constraint on the submatrix exponent value 
f(i,j) = g r {i) ■ ffc(j) (since all entries in each of g r and g c 
are distinct, |[T5l ). 



The SCB structure imposes certain conditions 1 15 1 on the 
variable and check nodes: 

Bit Consistency: The neighboring check nodes of a variable 
node must have distinct row-group (i) labels. 

Check Consistency: The neighboring variable nodes of a 
check node must have distinct column-group (j) labels. 

Pattern Consistency: (As shown in [15|) Since every entry 
in a row (i, I) and a column (j, k) in the SCB mother matrix 
with the value 1 satisfies k + ij = I mod p, if two variable 
nodes corresponding to columns (ji,kx) and (32^2) share a 
check node in row group i, they must satisfy: 



ki + iji — k 2 + ij2 mod p. 



(1) 



Converse of Pattern Consistency: If Q is satisfied, then the 
two variable nodes (ji, k\) and (J2, k 2 ) share a check node in 
row group i of the SCB mother matrix. 

Cycle Consistency: As shown in 1T51 . the equations of the 
form ([T} for any length-2< cycle in an SCB mother matrix, 
which involve t variable nodes with column-group labels j\ 
through jt and t check nodes with row-group labels i\ through 
it, show that the cycle must satisfy: 



E 



m \J (m+1 



mod t -jm) = mod p. 



(2) 



Converse of Cycle Consistency: If (|2]l is satisfied, a series 
of equations of the form ([TJ must also hold for some set 
{ki, . . . , kf} with < k m < p. For example, selecting any 
value of < k\ < p forces the values of {fc 2 , ■ ■ • , k t }. These 
equations of the form ([T| force the presence of the cycle by 
the converse of the pattern consistency condition above. 

One particular type of SCB code considered in this paper is 
formed by selecting a set of r rows from the SCB mother 
matrix. We call such codes selected-row (SR) SCB codes. 
Using the converse of cycle consistency, we can establish that 
all SR-SCB codes with r > 3 have a girth of 6 (i.e., the 
smallest cycles are always length-6 cycles.) 

Lemma 1: All SR-SCB codes with r > 3 have girth 6. 
Proof: It is known that the girth of such codes must be at least 
6. See O, 03). That the girth is exactly 6 follows from the 
fact that if r > 3, there is always a solution to 

ii0*2 - h) + iiih + htii ~ 33) = mod p , (3) 

which is ji = is, ji = H> an d J3 = ii- B 
Note that it may be possible to remove these length-6 
cycles by removing columns from the SR-SCB code (i.e by 
shortening the code). Although already established in [6|, 1 15 1, 
cycle consistency shows that length-4 cycles are not possible 
in any SCB codes because 

ii(j2 -3i) + 12(31 -32) = (k ~ i-2)(h - Ji) mod p. 

After reviewing absorbing sets below, Q is used to con- 
struct a matrix equation based on the cycles contained in an 
absorbing set that must be satisfied if that absorbing set is 
present in the SCB code. 
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Fig. 1. Depiction of a (4, 8) absorbing set. Black circles are the four variable 
nodes (bit nodes) of the absorbing set. The white squares are the satisfied 
neighboring check nodes, and the black squares are the eight unsatisfied 
neighboring check nodes. 

B. Applying the cycle consistency matrix to absorbing sets 

This section introduces the cycle consistency matrix as a 
way to capture the cycle behavior of a graphical structure. 
While the cycle consistency matrix approach is generally 
applicable to graphical structures that contain cycles in SCB 
codes, we provide the specific example of absorbing sets to 
illustrate the technique since the absorbing sets are the focus 
of the analysis in the rest of the paper. 

An LDPC code with parity-check matrix H is often viewed 
as a bipartite (Tanner) graph Gh = (V,F,E), where the set V 
represents the variable nodes, the set F represents the check 
nodes, and the set E corresponds to the edges between variable 
and check nodes. 

For a variable node subset V as c V, analogous to G#, let 
G as = (Xs.Fasj-Eas) be the bipartite graph of the edges E as 
between the variable nodes V as and their neighboring check 
nodes F- ds . Let o(V as ) C F as be the neighbors of V as with 
odd degree in G as and e(V. is ) C F w be the neighbors of V as 
with even degree in G as . We refer to the nodes in e{V a& ) as 
"satisfied check nodes" because they will satisfy their parity- 
check equations when all the nodes in are in error. 

Definition 2 (Absorbing Set (cf. H15V )): An (n v ,n ) ab- 
sorbing set is a set V RS c V with | V^ s | = n v and |o(T4s)| = n o» 
where each node in V as has strictly fewer neighbors in o(V as ) 
than in e(T4 s ). ■ 

Moreover, if each variable node in V\V as has strictly fewer 
neighbors in o(V as ) than in F\o(V as ), an (n v ,n ) absorbing 
set is called an (n v , n ) fully absorbing set [15 |. An important 
property of fully absorbing sets is that they are stable under 
bit-flipping decoding in which the bit values of variable nodes 
are flipped if a majority of their neighboring check nodes are 
not satisfied. 

Fig. [T] shows an example of a (4, 8) absorbing set, which 
has 4 variable nodes and 8 unsatisfied check nodes. 

Suppose there are n v variable nodes in an absorbing set. 
Let ji , . . . , j nv be the column-group labels of these n v nodes 
in the SCB mother matrix. Define 

u m = j m - ji for to e {2, n v } (4) 

and define u = [u2, ...,u riv }. For each cycle in the absorbing 



set, by replacing the difference of fs with the difference of 
it's, |2) may be written as 

t 

^ (im-i - im)u m = mod p, (5) 

m=2 

where 2t is the length of that cycle. Note that the sequence 
of check-node row groups {ii,i2, ■ ■ ■ , it} will be different for 
different cycles reflecting the particular cycle trajectories. 

Every cycle in the absorbing set satisfies an equation of the 
form (|5). Taken together, these equations produce a matrix 
equation: Mu = mod p, where M^ m is the coefficient of 
u m in |5]) for the £th cycle. 

Definition 3 (Cycle Consistency Matrix): A Cycle Consis- 
tency Matrix (CCM) M of a graphical structure containing 
cycles in a bipartite graph of an SCB LDPC code, for example 
an absorbing-set graph G as , has a number of columns one 
less than the number of variable nodes in the structure. Each 
row in M corresponds to a cycle in the structure, specifying 
the coefficients of u in (|5j for the corresponding cycle. 
Furthermore, Mu = mod p completely characterizes the 
requirement that every cycle in the structure satisfies ((5). ■ 

C. Efficient cycle consistency matrices 

A key property of a CCM is that Mu = mod p 
completely characterizes the requirement that every cycle in 
the structure satisfies ((5J. Even so, it is not necessary for M 
to include a row for every cycle in the structure. A cycle need 
not be included in M if it is a linear combination of cycles 
already included in M. Thus the number of rows needed in M 
is the number of linearly independent cycles in the structure. 

The notion from graph theory QUI of an incidence matrix 
is useful for establishing the number of linearly independent 
cycles in a graphical structure. 

Definition 4 (Incidence Matrix): For a graph with n ver- 
tices and q edges, the (unoriented) incidence matrix is an n x q 
matrix B with By = 1 if vertex Vi and edge Xj are incident 
and otherwise. ■ 

Note that since each edge is incident to exactly two vertices, 
each column of B has exactly two ones. 

The incidence matrix of a graph is useful for identifying 
the cycles in the graph because every cycle has the property 
that the indicator vector x c of the edges in the cycle satisfies 
Bx c = mod 2. In fact, the edges identified by the vector 
x c form a cycle (or a union of cycles) if and only if Bx c = 
mod 2. This is formalized in the definition below. 

Definition 5 (Binary Cycle Space): The binary cycle space 
(bcs) of a graph is the set of binary indicator vectors of the 
edges in a cycle or a union of cycles. The bcs is the null space 
over GF(2) of the graph's incidence matrix. ■ 

An efficient CCM is one with the least possible number of 
rows, which is the dimension of the binary cycle space £>bcs- 

Definition 6 (Efficient Cycle Consistency Matrix): An Effi- 
cient Cycle Consistency Matrix M of a graph G has £>bcs 
rows. The rows of M correspond to Db cs linearly independent 
cycles in G. Each row has the coefficients of u in |5]l for the 
corresponding cycle. ■ 
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For any bipartite graph, the check nodes and variable nodes 
together can be the nodes used to construct an incidence matrix 
whose nullspace is the binary cycle space. As an example, 
for the absorbing set of Fig. [JJ an 18 x 20 incidence matrix 
•8(4,8) can be formed from the 18 nodes (14 check nodes and 
4 variable nodes) and 20 edges in that graphical structure. The 
rank of -Bug) is 17 which indicates that D\, cs = 20 — 17 = 3. 
Thus there are three linearly independent cycles in the (4, 8) 
absorbing set of Fig.[T]and an efficient CCM for this absorbing 
set will have three rows. 



D. The VN graph approach to finding Df, cs 

The smallest absorbing sets for the SCB codes we study 
below all share the property that within the absorbing set, the 
check nodes all have degrees of two or less. This is consistent 
with work such as QTI in which elementary trapping sets, 
defined to have check node degrees of two or less l32l . are 
shown to dominate the code performance. 

The constraint that the check node degree be two or 
less permits a simpler technique for identifying the smallest 
number of rows needed for a CCM than what was presented 
above. Our examples of this technique are all absorbing sets, 
but the technique applies to elementary trapping sets and more 
generally any graphical structure in which the check nodes all 
have degrees of two or less. 

For any bipartite graph in which check nodes all have degree 
two or less a "variable-node" (VN) graph can be constructed 
whose only vertices are the variable nodes of the original graph 
and for which two vertices in the VN graph are connected 
iff there is a check node that connects them in the original 
graph. Under this constraint, for codes with girth greater than 
4 (including all SCB codes), multiple edges are not allowed 
between vertices of the VN graph since they would correspond 
to a length-4 cycle. 

As long as the check node degree is constrained to be two 
or less, the null space of the incidence matrix for the VN graph 
identifies the same binary cycle space for the absorbing set as 
the null space of the incidence matrix for the original bipartite 
graph. However, the VN graph greatly simplifies the process. 

The (4, 8) absorbing set provides an example of how to 
use the VN-graph approach to construct the CCM for a given 
absorbing set. The absorbing set graph in Fig. [JJ induces the 
variable-node (VN) graph shown in Fig. [2] There are five 
cycles in the variable-node graph, but not all of these cycles 
need to be explicitly represented in the CCM. The incidence 
matrix _B as of the VN graph shown in Fig. [2] with q — 6 edges 
and Vas — 4 nodes is 



1110 
10 110 
10 10 1 
10 11 



(6) 



The VN-graph incidence matrix _B as is simply the transpose 
of the submatrix H as of the parity-check matrix whose rows 
correspond to satisfied check nodes of G as (a subset of F as ) 
and whose columns correspond to the variable nodes in V as : 




Fig. 2. Variable-node (VN) graph of the (4, 8) absorbing set of Fig. [TJ 



110 
10 10 
10 1 



1 1 
1 1 



11 



(7) 



The rank of £? as in |6]i is 3 under GF(2). Thus the dimension 
of the binary cycle space is £>bcs = q— rank(£> as ) = 6—3 = 3, 
which means that three linearly independent cycles span the 
binary cycle space. This is the same answer we found above, 
but the VN graph approach uses a 4 x 6 incidence matrix rather 
than an 18 x 20 incidence matrix. 

Applying |5]) to the following three linearly independent cy- 
cles in Fig. [2] v\— v%— v^— V\, v\—V2—v&—V\, v%— V3— V4— v\ 
(here and in the remainder vi corresponds to the variable node 
labeled (jj, fcj) in Fig. [TJ produces the three rows of the CCM 
in ([8]i and its determinant in |9]) as follows: 





M = 



ll -l 2 l 2 - 15 

ii—i 6 i 6 - i 4 
«5 - h h~ U 



(8) 



and 
IMI 



-(ii-taXie-uX^— is) - (is— isXii— ieXis-u)- (9) 



E. Necessary and sufficient conditions for absorbing sets 

Recall that the elements of u as defined in |4]) contain 
difference information about the column groups involved in the 
absorbing set: the value of the first column group is subtracted 
from each of the others. The vector u cannot be an all-zero 
vector because the Check Consistency condition requires that 
variable nodes sharing a check node have distinct column 
groups, and a zero entry indicates that the variable node is 
in the first (j = ji) column group. Thus a necessary condition 
for the existence of a given absorbing set in an SCB code is 
that its M does not have full column-rank in GF(p). 

Definition 7 (Extensible VN Graph): If the VN graph of an 
absorbing set G as is a sub-graph of the VN graph of at least 
one other absorbing set G as with the same number of variable 
nodes, then the VN graph of G as is extensible. ■ 
Note that the incidence matrix B 2 of absorbing set G as is 
then the incidence matrix B\ of absorbing set G as with some 
additional columns. 

Equipped with these definitions, the next theorem gives nec- 
essary and sufficient conditions for the existence of absorbing 
sets in SCB codes. 
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Theorem 1: Given a proposed (n v , n Q ) absorbing set graph 
G as = 04sj -Pas, Eos), where every variable node is involved in 
at least one cycler] specified column group labels of the n v 
variable nodes in V as in the SCB mother matrix, and specified 
row-group labels of the check nodes in F as in the SCB mother 
matrix, the following are necessary conditions for the proposed 
absorbing set to exist in each daughter SCB LDPC code (with 
a parity check matrix H that includes the specified row and 
column groups of that SCB mother matrix): (1) The CCM 
for G as does not have full column-rank; (2) Variable nodes 
in Vas satisfy the Check Consistency condition and can form 
a difference vector u in the null space of the CCM; and (3) 
Each check node in F as satisfies the Bit Consistency condition. 
Taken together, these conditions are also sufficient if the VN 
graph of this absorbing set is not extensible. 

Proof: Each of the three conditions has already been shown 
above to be a necessary condition for the existence of G as in 
an SCB. Consider the sufficiency of the three conditions. 

If all three conditions are satisfied, all the cycles correspond- 
ing to the rows of the CCM exist in Gh because they can 
be constructed as follows: Conditions (1) and (2) ensure a 
sequence of variable node column groups \ji,j2, ■■■,3n v ] that 
form a vector [w 2 , ■■■,u nv ] m me null space of the CCM. For 
any fixed ki, we can compute k2,-~,k Uv using the converse 
of Cycle Consistency. Any linear combination of these cycles 
exists in Gh as well. 

These cycles cover every edge except edges that connect a 
variable node to a degree- 1 check node in G as . If (1) a variable 
node's unsatisfied check node is actually the same as another 
variable node's unsatisfied check node, or (2) a variable node's 
unsatisfied check node is actually the same as a satisfied check 
node in the graph, or (3) two of the satisfied check nodes are 
actually the same, then there exist other independent cycles in 
the VN graph which extend the VN graph. 

In these cases, the original structure is extensible and the 
three conditions are not sufficient to establish the existence 
of the absorbing set G as originally considered. The three 
conditions might instead be caused by the presence of an 
absorbing set whose VN graph contains the VN graph of the 
originally considered absorbing set as a subgraph. However, if 
the original VN graph is not extensible, the above constructed 
solution must establish the existence of the proposed absorbing 
set G as m Gjj. In this case the three conditions are sufficient. 
■ 

Corollary 1: Suppose a non-extensible (n Vl ,n 0l ) 
absorbing-set graph G as = (V^,F^,E^ S ) has a VN graph 
that is a sub-graph of the VN graph of another (n V2 , n Q2 ) 
absorbing-set graph G as = (V£, F a 2 s , E^). (Note that since the 
VN graph of G as is non-extensible, n Vl ^ n V2 .) Then, the 
existence of G as is a necessary condition for the existence of 

Gas- 

Proof: Suppose the CCMs of G^ and G^ are M 1 and M 2 
respectively. If the VN graph of G as is a sub-graph of the 
VN graph of G as , the independent cycles of G as will also be 
independent in G as and thus at least one valid M2 has Mi as 

2 If the variable node degree is at least 2, then each variable node in a given 
absorbing set must be a part of at least one cycle. 



a sub-matrix: 



Mo 



Mi 
A B 



(10) 



where the sub-matrix [A B] represents the other linearly 
independent cycles in G as , which are not included in G as . 
Therefore if there exists a valid 112 such that M2U2 = 
mod p, the first n v \ — 1 elements would also be a valid Ui 
such that M1U1 = mod p. Because G as is non-extensible, 
a valid Ui such that M1U1 = mod p is sufficient to force 
the existence of G as . This shows that the existence of G as is 
a necessary condition of the existence of G as . ■ 



IV. Analytical results for r = 5 SCB codes 

We now consider the H^fu -\ LDPC codes as an instance 
of high-performance practical LDPC codes with a hardware- 
friendly structure. We consider r = 3, r = 4 and r = 5 SCB 
codes, but we begin with the r = 5 case because simple row 
selection eliminates the smallest absorbing sets in this case. 
Thus, this section considers r = 5 (five row groups) and shows 
how to design an SCB code with a specified circulant matrix 
that eliminates the smallest absorbing sets by selecting rows 
from the SCB mother matrix to force the CCMs associated 
with the dominant absorbing sets to be full column-rank. 

This example of SCB code design involves two classes of 
SCB codes: 

• Array-based codes |6) are the most elementary SCB 
codes in which the first r rows of the SCB mother matrix 
Hp' P f : f{i,j) = i • j comprise the parity-check matrix. 
We will refer to this class as the elementary array-based 
(EAB) codes. 

• As shown in ifTOl . a careful selection of the r row- 
groups from the overall SCB mother matrix can improve 
performance over the EAB codes. Thus, selected-row 
(SR) SCB codes are our second class of SCB codes. 
The parity-check matrix for these codes is H^, f(i, j) = 
g r (i) ■ j where g r (i) is called the row-selection function 
(RSF). We will often represent an RSF as the vector 
MO) 3,(1) g r (2) g r {3) <? r (4)]. 

Theorem [T] shows that an absorbing set may be avoided 
either by forcing the associated CCM to have a zero-dimension 
nullspace or by precluding u from being in the null space 
of M. Corollary [T] shows that if a non-extensible absorbing 
set does not exist, then all absorbing sets whose VN graphs 
contain the VN graph of this absorbing set also do not 
exist. The CCM approach carefully selects the RSF (and 
possibly also an analogous column selection function (CSF)) 
to preclude small absorbing sets, in the order of the size (the 
value of n v ) of the VN graph of the absorbing sets. 



Section |IV-A establishes that (4, 8) absorbing sets are the 



smallest possible for a general r = 5 code family. Then 
Corollary [2] shows that (4, 8) absorbing sets indeed exist for 
this r = 5 array-based code family for each set of row 
group labels that satisfy the Bit Consistency conditions and 
dct M — mod p. This theoretical result is also consistent 
with previous experimental results of a sum-product decoding 
algorithm implemented in software and on a hardware emu- 
lator |3J for which it was shown that decoding errors due to 
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(4, 8) absorbing sets dominate the low BER region of certain 
r = 5 array-based codes. 

Sections |IV-A| |IV-B and IV-C respectively, show that 



(4, 8), (5, 9) and (6, 8) absorbing sets exist in the EAB code 
with r = 5 and that carefully selecting row groups from 
the SCB mother matrix can preclude all these absorbing 



sets. Section IV-D provides several good RSFs that preclude 
these absorbing sets for r — 5 SCB codes. Section |IV-E| 
explores the absorbing set spectrum of the existing quasi-cyclic 
LDPC codes with the Tanner construction [4|. Section [IV-F| 
identifies equivalence classes among SCB codes and uses the 
equivalence-class approach to enumerate all of the RSFs for 
r = 5 and p = 67 that preclude all (4, 8), (5, 9) and (6, 8) 
absorbing sets. 

A. (4, 8) absorbing sets 

In this section we analyze (4, 8) absorbing sets. The main 
result states that (4, 8) absorbing sets exist in SCB codes 
and specifically in EAB codes (Corollary [3j, but that (4, 8) 
absorbing sets can be provably eliminated from SR-SCB 
codes using a suitable RSF, as shown in Corollary |4] These 
theoretical results will be substantiated by experimental results 
in Section IVIII 

Since Lemma[T]establishes that all SR-SCB codes have girth 
6, (4, 8) absorbing sets are the smallest possible for SR-SCB 
codes with variable node degree r = 5. (Also see |9).) To 
see why this is true, note that the cases with n v < 4 are 
not possible because each of n v variable nodes needs to have 
at least three satisfied checks but this would require a girth 
of four. When n v = 4, the girth constraint and the need 
for three satisfied checks prevent check node degrees larger 
than two and prevent any variable node from sharing more 
than three checks with other variable nodes in the absorbing 
set. Hence, n — 8 and each pair of variable nodes in the 
absorbing set shares a distinct satisfied check. Therefore, the 
(4, 8) absorbing set in Fig.[T|is the smallest possible absorbing 
set in the H ,p fr > code family. More general results regarding 
the minimality of absorbing sets are provided in ||9J. 

This section shows that the (4, 8) absorbing set always exists 
in r = 5 EAB codes. However, by carefully selecting row 
groups from the SCB mother matrix, r = 5 SR-SCB codes 
can preclude all (4, 8) absorbing sets. 

Using Fig. [T] and |9]), Theorem [T] leads to the following 
corollary regarding (4, 8) absorbing sets: 

Corollary 2: Consider an SCB mother matrix with a spec- 
ified p. The existence of a selection of integers for row group 
labels ii,. . . ,ig satisfying the Bit Consistency conditions as- 
sociated with the absorbing set shown in Fig. fT] and satisfying 
det M = mod p (with det M given in is necessary 
and sufficient for the existence of the (4, 8) absorbing set of 
Fig. [T]in the bipartite graph of the SCB mother matrix. 

Proof: Since the VN graph in Fig. [2] of the (4, 8) absorbing 
sets is a fully connected graph, it is not extensible without 
introducing a parallel edge which implies a length-4 cycle in 
the corresponding bipartite graph. This property would violate 
Lemma[T] Because the VN graph is non-extensible, Theorem[T] 
applies. Hence dot M=0 mod p and Bit Consistency are both 



necessary conditions. If the Check Consistency is also satis- 
fied, the three conditions together are sufficient by Theorem 
1. The rest of the proof shows that identifying a selection 
of integers for row group labels ix, satisfying the 

Bit Consistency conditions and satisfying det M=0 mod p 
implies the existence of jx, 32,33, ji satisfying the Check 
Consistency conditions. 

Consider the Bit Consistency conditions associated with 
each of the four variable nodes in Fig. [T] Our concern is 
only regarding Bit Consistency conditions that involve satisfied 
check nodes. Hence the six row-group labels of interest are 
{ii,..., ig}. For example, the Bit Consistency conditions 
applied to (jx,ki) require that ix^ii, and £47^5. 

There are six different row-group labels of interest, but only 
five possible row groups since r = 5. At least one pair of row- 
group labels must share the same row group. The inequalities 
implied by the Bit Consistency conditions allow only i\ = 13, 
i% = ?4, and 15 = Iq. 

As shown in lfl5l for (4, 4) absorbing sets, if all three of 
these equalities are satisfied, the Cycle Consistency condi- 
tions applied to the three length-8 cycles involving two pairs 
of same-row-group check nodes (i.e., v\— V2— V3— U4— Vx, 
Vx—V2—V4,—V3—Vx, and v\— v±— v^— U3— v±), yield three 
equations that require all variable nodes to have the same 
column group. This violates the Check Consistency conditions. 

Visualizing the (4, 8) absorbing set as a triangular-base 
pyramid with the four variable nodes as vertices, every point 
is symmetric to every other point and the three possible 
equalities allowed by the Bit Consistency conditions are all 
isomorphisms. Thus, we only need to consider one case of 
one equality being satisfied and one case of two equalities 
being satisfied. 

Thus, there are only two possible non-isomorphic row- 
group labelings for the (4, 8) absorbing set as follows: 
for any bijective assignment of five distinct row group 
labels to {t,x,y,z,w}, (ix, i%, 13, U, 15, i&) can be either 
(x, y, x, y, z, w) (Assignment 1, where two equalities are satis- 
fied) or (x, t, w, y, z, z) (Assignment 2, where only one equal- 
ity is satisfied). A detailed proof of the possible assignments 
can be found in ifTOl . Applying (|9) to these two assignments 
yields the following: 

detM = (y - x) ((z - x)(w - y) + (z - y)(w - x)) (11) 

for Assignment 1, and 

detM = (z-w)(x-t)(y-z)-(y-w)(x-z)(z-t) (12) 

for Assignment 2. 

Suppose det M = mod p. Then there exists a non-zero 
solution to M • u = mod p, where u=[u2, U3, «4] T - With 
either Assignment 1 or Assignment 2, the resulting M has six 
nonzero entries such that if one of {u2,u 3 ,u 4 } is nonzero, 
all three must be nonzero. Furthermore, the structure of M 
ensures that any nontrivial solution to M • u = mod p 
will have u%, U3, and U4 all distinct. Thus, for a fixed jx, 
we can find j'2, J3, and j'4 without contradiction to Check 
Consistency (since all are distinct). Fixing one specific k 
value will determine all others according to Cycle Consistency, 
yielding the variable nodes of (4, 8) absorbing sets in the code. 
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Therefore, a row-group labeling that satisfies Bit Consistency 
and det M=0 mod p is a necessary and sufficient condition 
for the existence of (4, 8) absorbing sets. ■ 

We will examine how Corollary [2] applies to all 
SR-SCB codes. To begin, consider the special case 
of EAB codes, which select the first five rows so 
that {x,y,z,w,t} = {0,1,2,3,4}. For Assignment 1 
there are no solutions that achieve det M=0 mod p 
for prime p large enough (p>17). However, for 
Assignment 2 there are 8 solution sets for (x,y,z,w,t): 
{(4, 3, 2, 0, 1), (4, 1, 2, 0, 3), (3, 4, 2, 1, 0), (3, 0, 2, 1, 4), 
(1, 4, 2, 3, 0), (1, 0, 2, 3, 4), (0, 3, 2, 4, 1), (0, 1, 2, 4, 3)}. Note 
that in these solutions z, the value of i$ — i§, is always 2. 

Remark 1: The eight Assignment 2 solutions above imply 
a relatively simple structure. The six edges of the VN graph 
in Fig. [2] include three pairs of edges such that the edges in 
each pair do not share any variable nodes. Each pair of edges 
corresponds to a pair of check nodes in Fig. [2] The check-node 
pairs are (i\, 13), (i 2 , h), and (is, is)- Including isomorphisms, 
the eight solutions above imply that one of the pairs of check 
nodes has both check nodes in row group 2, a second pair 
of labels must have check nodes in row groups and 4, and 
the third pair must have check nodes in row groups 1 and 
3. The eight solutions above assume that the (k,k) pair has 
both check nodes in row group 2, but choosing one of the 
other pairs to have both check nodes in row group 2 merely 
produces an isomorphism. 

As an example, let's find one set of j's and fc's so that 
det M=0 mod p. Among the EAB Assignment-2 solutions 
above, consider the last one: (x=0, y=l, z=2, w=4, t—3) so 
that (ii,k,k,k,k,k) = (x,t,w,y,z,z) = (0,3,4,1,2,2). 
For this solution set, 
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and 



detM = _(-3)(l)(-2) - (l)(-2)(3) = . (14) 

As shown in ( fT4j >, det M=0 for any p. Suppose p = 47 
and u 2 — 1. Solving Mu = gives us u 3 — 3 and U4 = 2. 
Selecting j\ = gives j 2 = 1, j'3 = 3,j4 = 2. From Pattern 
Consistency, taking fci = gives k 2 — 0, k 3 — 41,^4 = 
45. Thus, once a nonzero solution to Mu = mod p is 
specified (for example by selecting u 2 and solving for u as 
above), specifying any ji and ki leads directly to a solution 
that identifies a specific absorbing set. 

Continuing with this specific choice of (ii, i 2 , k, k, k, k)< 
there are p—1 ways to select u 2 , each of which yields a distinct 
nonzero solution to dctM=0 mod p. (Selecting u 2 =0 gives 
the all-zeros solution.) For each of these solutions there are p 
ways to select each of ji, and k%, yielding a total of p 2 (p— 1) 
solutions. 

Thus there are there 0(p 3 ) such solutions. Here and in 
the remainder of the paper we use the standard asymptotic 
notation for 0: a positive function f(p) is &(p k ) if there exist 
constants c\ and c 2 such that < c\ < c 2 < 00 for which 



CiP k < f(p) < c 2 p k for all p > 0. For each of the eight 
Assignment-2 solutions above there are ®(p 3 ) absorbing sets. 
The overall number of absorbing sets is thus 0(p 3 ). 

To avoid the (4, 8) absorbing set we need to consider 
alternatives to the EAB selection of row-group labels in 
{0,1,2,3,4}. In general, SR-SCB codes can select the row- 
group labels to force the CCM for the (4, 8) absorbing sets 
to have a zero-dimension nullspace (equivalent to forcing the 
determinant to be nonzero for this square CCM). The goal 
is to identify a set of five (since r = 5) row-group labels 
such that any valid assignment of those labels to i\ through 
ig satisfies det M^0 mod p in both ( fTT) and ( |12) . One such 
example is the RSF [0,1,2,4,6]. For a valid assignment of 
{0,1,2,4,6} to i\ through k, recall that the Bit Consistency 
conditions allow only i\ — 13, i 2 — k, and/or i$ — k- Any 
such valid assignment of these labels has det M^0 mod p in 
both (jTTJ and ( fT2| ) for any prime p greater than 23. Therefore 
in SR-SCB codes with a well chosen RSF there are no (4, 8) 
absorbing sets for sufficiently large p. Another example is the 
RSF [0,1,3,8,19] which avoids all (4,8) absorbing sets for 
p — 47 and certain other p's as described in iflOl . 

The above analysis proves the following corollaries. 

Corollary 3: The (4, 8) absorbing sets exist in all r = 5 
SCB mother matrices and in particular in all EAB codes de- 
scribed by the parity check matrix H and their multiplicity 
scales as 0(p 3 ) in the EAB codes. 

Corollary 4: There are no (4, 8) absorbing sets in the SR- 
SCB codes described by the parity check matrix H ' p ,., ., for 
prime p large enough and with a proper choice of g r (i). 

B. (5,9) absorbing sets 

Assuming an RSF that avoids the (4, 8) absorbing sets, this 
section proves that the (5, 9) absorbing sets are the smallest 
remaining. The CCM approach shows that the (5, 9) absorbing 
sets always exist in the EAB codes (Corollary [8j, but are 
avoided for SR-SCB codes by some of the RSFs that precluded 
the (4, 8) absorbing sets (Corollary ]9J. We start by establishing 
a series of intermediate results, proven in Corollaries [5] [6] 
and |7l 

Corollary 5: For an r = 5 SR-SCB code in which (4, 8) 
absorbing sets do not exist, (5, b) absorbing sets do not exist 
for 6 < 9. 

Proof: For r = 5, the total number of edges in a (5, b) 
absorbing set is 25. The number of edges that connect to 
satisfied check nodes is even, with two such edges associated 
with a pair of connected variable nodes. Since the girth is 
six for SR-SCB codes, each pair of variable nodes can be 
connected at most once. There are ten such pairs. Any group 
of four variable nodes forms a (4,8) absorbing set iff all six of 
its variable-node pairs are connected. Thus, if (4, 8) absorbing 
sets do not exist, each such group of four variable nodes must 
have one unconnected pair. This implies that at least two of the 
ten possible pairs must be unconnected. Thus, there are at most 
eight connected variable-node pairs, which implies at most 
sixteen edges connected to satisfied check nodes, which means 
at least nine edges connected to unsatisfied check nodes, and 
therefore b > 9. ■ 
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Fig. 3. Depiction of a (5,9) absorbing set. Black circles are bit nodes in 
the absorbing set, white squares are their satisfied checks, and black squares 
are their unsatisfied checks. 

Corollary 6: For an r — 5 SR-SCB code in which (4, 8) 
absorbing sets do not exist, the (5, 9) absorbing set is the 
smallest possible, and there are no other (5,6) absorbing sets. 

Proof: Because each of the five variable nodes must have 
three satisfied checks, fifteen of the twenty five edges must 
go to satisfied checks. Thus, the number of edges that go to 
unsatisfied check nodes is at most ten. From Corollary [5] 6>9 
in the absence of (4, 8) absorbing sets. Since 6 is odd, the 
smallest possible (5,6) absorbing set is the (5,9) absorbing 
set. We now prove that the (5, 9) absorbing set does not contain 
the (4, 8) absorbing-set as a sub-graph. 

Suppose a (5,9) absorbing set has a (4,8) absorbing set as 
a subgraph. We label the five variable nodes as Ui, U2j V3, W4 
and v$. Without loss of generality, let vi,V2,vs and V4 form 
a (4, 8) absorbing set as in Fig. [TJ If v$ connects to a satisfied 
check node in the (4, 8) absorbing set, for example check node 
i\ that connects v\ and v%, v% and v 2 will have more unsatisfied 
check nodes than satisfied check nodes. Then 175 has also to 
connect to one of vi's unsatisfied check nodes (i-j or in) to 
make v\ satisfy the absorbing set property. However, Vx,ii, v& 
and 17 or Vi,i%, v$ and in would then form a length-4 cycle, 
which leads to a contradiction of Lemma [TJ 

Thus V5 cannot connect to the satisfied check nodes of 
v i,V2,vs and V4, Since t> 5 has at least three satisfied check 
nodes, it has to connect to at least three of the formerly 
unsatisfied check nodes of V\,V2,v§ and V4. This makes 
the total number of unsatisfied check nodes seven or less 
and certainly not nine. This again leads to a contradiction. 
Therefore (5, 9) absorbing sets cannot have (4, 8) absorbing 
sets as a subset. ■ 

The only remaining possible configuration of a (5, 9) ab- 
sorbing set is where one variable node has four satisfied check 
nodes and four variable nodes have three satisfied check nodes. 
Fig. [3] depicts this configuration, which has a VN graph that 
does not contain the VN graph of a (4, 8) absorbing as a 
sub-graph. Corollary [7] below establishes the necessary and 
sufficient conditions of this absorbing set. 

Corollary 7: For matrix M in ( p"5j ) below, a row-group label 
assignment that satisfies the Bit Consistency conditions and 
satisfies det M = mod p is necessary and sufficient for 
the existence in an SR-SCM code of the (5,9) absorbing set 



shown in Fig. [3] 

Proof: With an analysis similar to that of (4, 8) absorbing 
sets, the binary cycle space for Fig. [3] has dimension 4. 
We construct the following CCM by selecting the following 
linearly independent cycles: 
V1—V3—V4—V1, and V\— U4- 



M = 



1 1 
h 



«5 



»5 - 22 



*2 - i& 




v 1 -v 2 -v 3 -v 1 , v 1 -v 2 -v 5 -v 1 , 

V 5 -Vi. 




H - h 
h - h 



i 8 - 14 




i 7 - ^4 



(15) 



As in the proof of Corollary [2] we can show that a 
row-group label assignment that satisfies the Bit Consistency 
conditions and dct M = mod p is necessary and sufficient 
for the existence of a (5,9) absorbing set, where 



detM =(ii - i 5 )(i 8 - i 4 )(i 2 - ieX^ - h) 

- (k - «8)(«5 - «2)(«6 - is)(*7 - U) 



mod p. 
(16) 



Under the Bit Consistency conditions, there are five 
possible non-isomorphic check labeling patterns of 
(ii,i2,*3)*4,«5)*6i*7j*8) as follows: (x,y, z,w, z,w,x,y), 
(x, y, z, w, t, w, x, y),(x, y, z, w, t, w, x, z),(x, y, z, w, t, x, t, y), 
and (x,y,z,w,t,w,t,y), where different letters correspond 
to distinct row -group labels. 

The EAB codes have the set of row- 
group labels {x, y, z, w, t} drawn from the set 
{0,1,2,3,4}. For the 4th pattern (x,y,z,w,t,x,t,y), 
there are eight solution sets (x, y, z, w, t) e 
{(4, 0, 1, 3, 2), (4, 0, 3, 1, 2), (3, 1, 4, 0, 2), (3, 1, 0, 4, 2), 
(1, 3, 4, 0, 2), (1, 3, 0, 4, 2), (0, 4, 1, 3, 2), (0, 4, 3, 1, 2)} that 
always have det M=0 in ( [To} . The other four patterns have 
a nonzero determinant for p large enough. Once the labels 
of the check nodes are selected so that det M=0 in ( fTo} 
(cf. Fig. [5]), the variable node labels (pairs (ji,ki) through 
(j5,/c 5 )) can be selected in 0(p 3 ) ways as in the (4,8) case. 

For SR-SCB codes with a proper choice of RSF, det M^0 
mod p for p large enough. One such example is using an 
RSF of [0,1,2,4,7] where detM^O mod p for the prime 
p greater than 89 except for the set p e {101, 103, 131, 179}. 
Therefore we can conclude with the following two corollaries. 

Corollary 8: (5, 9) absorbing sets exist in EAB codes de- 
scribed by the parity check matrix H'?*, and their number 
scales as 0(p 3 ). 

Corollary 9: There are no (5,9) absorbing sets in the SR- 
SCB codes described by the parity check matrix H 5,p , ., ., for 
prime p large enough and with a proper choice of RSF g r (i). 

C. (6, 8) absorbing sets 

This section considers the (6, 8) absorbing sets, which are 
the smallest remaining after the (4,8) and (5,9) absorbing 
sets. We will investigate the (6, 8) absorbing sets both for 
EAB codes and for SR-SCB codes that preclude the (4,8) 
and (5, 9) absorbing sets. The six subsections that follow 
examine respectively the six candidate configurations of (6, 8) 
absorbing sets to be studied. 
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Section summary. We first prove in subsection |I V-C 1 1 that 
in (6,8) absorbing sets in r = 5 SR-SCB codes that preclude 
(4, 8) absorbing sets the largest possible check node degree 
is 2. Combinatorial and consistency arguments show (sub- 
sections IV-C3 and IV-C4i that three of the remaining five 



configurations are not present for p sufficiently large in either 
the EAB code or in SR-SCB codes that preclude the (4, 8) 
and (5,9) absorbing sets. 

The other two configurations have the cardinality 0(p 3 ) in 
the EAB code. However, both of these configurations contain 
a (4, 8) absorbing set as a subset and thus cannot be present 
in SR-SCB codes that preclude the (4, 8) absorbing sets. 



These two configurations are analyzed in subsections IV-C2 
and lIV-C5l 

With this comprehensive analysis of (6, 8) absorbing sets we 
can conclude that r = 5 SR-SCB codes can avoid all (4, 8), 
(5, 9), and (6, 8) absorbing sets for p large enough if the five 
row groups are properly chosen. 

1) Configuration 1 - Check nodes with degree>2: This case 
is precluded by the following lemma: 

Lemma 2: If a (6, 8) absorbing set has a check node that 
connects to more than two variable nodes in the absorbing set 
graph, it must contain a (4, 8) absorbing set as a sub-graph. 

Proof: Consider a check node that connects to more than 
two variable nodes in such an absorbing set. It is either a 
satisfied check connected to either 4 or 6 variable nodes or an 
unsatisfied check connected to either 3 or 5 variable nodes. 

If this check is satisfied, it cannot be connected to all 6 
variable nodes because any additional satisfied check would 
complete a length-4 cycle and violate the girth constraint. 
Suppose now that the check in question is satisfied and 
connected to 4 (out of 6) variable nodes. As with the degree-6 
case discussed above, these 4 variable nodes (forming a 4- 
clique) cannot have any other checks in common because of 
the girth constraint. 

The nodes in this 4-clique each have at least 3 satisfied 
checks, by the absorbing set constraint. For each of these 4 
variable nodes, two satisfied (degree-2) checks must necessar- 
ily be shared with the remaining 2 variable nodes. By the girth 
constraint, these 4 variable nodes then have exactly 3 satisfied 
and 2 unsatisfied checks each. The two remaining variable 
must share an additional satisfied check (so that all their checks 
are satisfied) to make the total number of unsatisfied checks be 
equal to 8. Consider two variable nodes from the 4-clique and 
the two remaining variable nodes outside of the clique. Each 
of these 4 variable nodes pairwise shares a satisfied check with 
the other 3 variable nodes. These checks are distinct and the 
configuration induced by these 4 variable nodes is indeed a 
(4, 8) absorbing set. 

Now, we suppose that the check connected to more than 
2 variable nodes is unsatisfied. It cannot be connected to 5 
variable nodes because to have three satisfied checks, pairs of 
these five variable nodes would need to share at least one ad- 
ditional check, violating the girth constraint. Suppose that this 
check is connected to 3 variable nodes. These 3 variable nodes 
(creating a 3-clique) must each share a distinct satisfied check 
with each of the remaining 3 variable nodes, by the absorbing 
set property. By the girth constraint, each variable node in the 




Fig. 4. Configuration 2 of a (6, 8) absorbing set. In this configuration two 
variable nodes have 5 satisfied check nodes which forces the remaining four 
variable nodes each to have exactly two satisfied check nodes. Note the two 
(4,8) absorbing sets (wi , V2 , V3 , v±) and (vz, U4, vs,v%), 

clique has an additional unsatisfied check not shared with any 
other variable node. The total number of unsatisfied checks 
incident to the 3-clique is thus 4. Since the total number of 
unsatisfied checks is 8, the remaining 3 variable nodes likewise 
share one unsatisfied check (thus creating another 3-clique) 
and each have an additional unsatisfied check not shared with 
any other variable node. Consider two variable nodes from 
each of the two 3-cliques. In this group of 4 variable nodes, 
every pair of variable nodes is connected via a distinct check, 
forming a (4, 8) absorbing set that is a subgraph of the (6, 8) 
absorbing set. 

Hence, all cases have been considered, and any (6, 8) 
absorbing set with a check node having degree larger than 
two must contain a (4, 8) absorbing set as a sub-graph. ■ 

Based on Lemma [2] attention is now restricted to cases 
where all check nodes in the absorbing set graph have degree 
at most 2. 

In a candidate (6, 8) absorbing set, variable nodes can have 
3, 4 or 5 satisfied checks. To maintain 8 unsatisfied checks, 
there can be at most 2 variable nodes with 5 satisfied checks. 
Suppose first that there are two such variable nodes. Since 
there are a total of 8 unsatisfied checks, the other 4 variable 
nodes must each have 3 satisfied and 2 unsatisfied checks. This 
necessarily implies the configuration shown in Fig. [4] which 
we discuss next. 

2) Configuration 2 - Fig. |5J The configuration in Fig. [4] 
has two overlapping (4, 8) absorbing sets. Thus the existence 
of this (6, 8) absorbing set requires zero determinants for 
the CCMs of these two (4, 8) absorbing sets. Thus, SR-SCB 
codes that preclude (4, 8) absorbing sets also preclude this 
(6, 8) absorbing set configuration. This is formalized by the 
following corollary. 

Corollary 10: There are no (6, 8) absorbing sets with the 
structure of Fig. [5] in SR-SCB codes that preclude (4, 8) 
absorbing sets. 
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We can also show the following result: 

Corollary 11: The existence of a selection of integers for 
row group labels satisfying the Bit Consistency conditions as- 
sociated with the absorbing set shown in Fig. [4] and satisfying 
det Mi = mod p and dct M2 = mod p, where Mi 
and M2 are CCMs of the two internal (4, 8) absorbing sets 
constitute a necessary and sufficient condition for the existence 
of the (6, 8) absorbing set in Fig. [4] 
The proof of Corollary [TT| follows that of Corollary [2] 

For EAB codes, we can show the following striking result: 

Corollary 12: For EAB codes, there are exactly p 2 (p — 1) 
(6, 8) absorbing sets of the type shown in Fig. [4] for p > 17. 

Proof: The (6, 8) absorbing set in Fig. |4] includes two (4, 8) 
absorbing sets, one spanning the 4 variable nodes (ji,ki), 
(h, k 2 ), C?3, h) and (j 4 , fc 4 ) and the other spanning (j 3 , fc 3 ), 
(34, h), {35, fo) and (j 6 ,k e ). 

Recall from Remark [I] in Section IV-A that any (4, 8) 
absorbing set in an EAB code requires two check nodes 
connecting two disjoint pairs of variable nodes to be both in 
row group 2. Another two check nodes connecting two disjoint 
pairs of variable nodes need to be in row groups 1 and 3. The 
third and final two check nodes connecting two disjoint pairs 
of variable nodes need to be in row groups and 4. 

Applying this remark simultaneously to the two (4, 8) 
absorbing sets in Fig. [4] reveals that the only possibility that 
does not violate the Bit Consistency condition is for i\, 13, and 
ii5 in Fig. [4] to be in row group 2. Neglecting isomorphisms, 
with i\ — iz — H5 = 2, the required structure permits only 
two distinct row group labelings for the remaining satisfied 
check nodes. Both have £4 = ig = 1 and i 2 = iio = 3. 
However, one labeling has is = 114 = 4 and ig = 213 = 
while the other has i$ = in = and i$ = iis = 4. 

For each of these two isomorphisms, there are p 2 (l — p)/2 
ways to distinctly select u 2 , ji, and k\ (Symmetry causes half 
of the p 2 (l— p) choices for u 2 , ji, and k\ to produce a labeling 
that is a 180 degree rotation of the labeling resulting from a 
different choice of u 2 , ji, and k\.) Thus the total number of 
absorbing sets of the type in Fig. [4] is exactly p 2 {p — 1). ■ 

Suppose now that there is exactly one variable node in the 
absorbing set having all five checks satisfied. The variable 
nodes in the absorbing set must necessarily be arranged either 
as in Fig. [5] or Fig. [6] 

3) (6,8) configuration 3 - Fig. [3J For this configuration, 
we have the following result. 

Lemma 3: In the Tanner graph corresponding to j\ 
for the EAB codes and SR-SCB codes, there are no (6, 8) 
absorbing sets for p large enough of the type shown in Fig. [5] 

Proof: Without loss of generality we may assign check node 
labels for the checks emanating from the variable node (ji, ki) 
as follows: i\ = x,i 2 = V, 13 = z, ia — w, and 25 = t, where 
x, y, z, w, t are the five distinct check labels. 

The binary cycle space for Fig.|5]has dimension 6. We con- 
struct the CCM of ( fT7| i by selecting the linearly independent 
cycles v%— v 2 — vq—v\, vi-v^-vq-vi, vi-v 2 -v$-vi, vi-v^-v^-vi, 
V1-V3-V5-V1, and vi-v 2 -V4-vi as follows: 




Fig. 5. Configuration 3 of a (6, 8) absorbing set. This is the first of two 
configurations with exactly one variable node that has 5 satisfied check nodes. 



M 



x — iq 

y-i 9 

x-i 7 

y -in in ~z 



i 6 -t 

i 9 - 1 
i*! — w 

y - i w i 10 - w 
x — is ig, — z 



(17) 



(18) 



The rank of the matrix is at most 5 so in fact we may 
consider the 5x5 submatrix formed by the first five rows 
(call it B). If the matrix B is full rank, then M is full rank 
and has a zero-dimension nullspace. Hence det(S) = is 
necessary for the existence of the absorbing sets of this type. 
This condition can be expressed as 

-(ill ~ z)[~{x - i 6 )(i 9 - t)(i 7 - w)(y - i w ) 
+{x - i 7 ){ie - t){y - i 9 )(iio —w)] = mod p. 

Also consider the bottom-left 4x4 submatrix (call it A). 
Note that the Bit Consistency conditions ensure that every 
element of M that is not explicitly zero in ( fT7j ) must be 
nonzero. Thus, if the matrix A is full rank, then M is full 
rank (rank-5) because either of the top two rows provides a 
row linearly independent from the bottom four rows. Hence 
dct (A) = is necessary for the existence of absorbing sets of 
the type shown in Fig. [5] This condition can be expressed as 

-(% - ir)(iio ~w){y- in) (is - z) 
+ (x - i 8 )(i 7 - to) (in - z)(y - i w ) = mod p. 

For the values of ig, 17, ig, iio and in in the available label set 
{x,y,z,w,t} both for the EAB and for the SR-SCB codes, 
and such that the i labels of check nodes sharing a variable 
node are distinct (see Fig. [5]), ([18) and ( fT9] > evaluate to zero 
for only a finite number of values of the parameter p. 

For g r (i) = i (the EAB code) det(B) ^ for p > 23. 
We can also find many SR-SCB codes where ( fT8| l and ( fl"9| ) 
evaluate to zero only for a finite number of values of the 
parameter p. For example, with the SR-SCB code defined by 



(19) 
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the RSF [0,1,2,4,7], det (A) ^ for p > 89. Thus for p 
sufficiently large, neither the EAB code nor well-designed SR- 
SCB codes contain (6, 8) absorbing sets of the type shown in 
Fig.0 ■ 

Incidentally, there do exist values of p for certain RSFs for 
which both det(A) = and det(B) = 0. For example, for 
p = 11, 17, 19, 23 with the RSF [0, 1, 2, 3, 4] both conditions 
hold. As we show in Corollary 13 below, det (A) = and 
det(-B) = together are a sufficient condition for the existence 
of absorbing sets of the type shown in Fig.|5]if (4, 8) absorbing 
sets have been excluded by the RSF. We also identified RSFs 
such as [0, 1, 2, 4, 8] where both determinants evaluate to zero 
(even before applying the mod p) for some choice of check 
node labels, and thus the associated codes have these (6, 8) 
absorbing sets regardless of the choice of p. 

Corollary 13: In the Tanner graph corresponding to 



-ff 5 '?/. _-\, if (4,8) absorbing sets have been excluded by the 



RSF, (6,8) absorbing sets of the type shown in Fig. [3] exist if 
and only if row labels can be specified to satisfy Bit 

and Check Consistency conditions and to satisfy det (A) = 
and det(B) = 0, where A and B are as defined in the proof 
of Lemma [3] 

Proof: We already established above that both det (A) — 
and det(i?) = are necessary for the existence of (6, 8) 
absorbing sets of the type shown in Fig. [5] We only need to 
establish that these two conditions are also sufficient. 

The rank of any matrix is lower bounded by the rank of any 
of its submatrices. Recalling that every element of M that is 
not explicitly zero in ( flT) must be nonzero, C is a full-rank 
4 by 4 matrix. Thus the rank of M is no less than 4. 

The first three rows of A are always linearly independent. 
(Recall that every element of M that is not explicitly zero in 
( fP7| i must be nonzero.) Thus if det(A) = the bottom row 
of A is a linear combination of the other 3 rows. Thus, the 
bottom row of M is a linear combination of the three rows of 
M just above it, and rank(B) — ranfc(M). 

Thus, if dct(yl) and det(-B) are zero, M is not full rank, 
and there is a non-zero u in the nullspace of the M. If the 
selected RSF precludes (4, 8) absorbing sets, then the (6, 8) 
absorbing set in Fig. [5] is not extensible. Thus by Theorem [T] 
det(A) = and det(£>) = is a sufficient condition for the 
existence of the (6, 8) absorbing set in Fig. [3] ■ 

4) (6, 8) configuration candidate 4 - Fig. [6j Note that the 
configuration shown in Fig. |6] contains a subgraph that is the 
(4, 8) absorbing set of Fig. fTT This implies the following: 

Corollary 14: For SR-SCB codes that do not contain (4, 8) 
absorbing sets, the configuration of Fig. [6] is not possible. 

Even though EAB codes contain (4, 8) absorbing sets, 
the following lemma states that EAB codes avoid the (6, 8) 
absorbing set in Fig. [6] 

Lemma 4: In the EAB codes corresponding to H 'f ., there 
are no (6, 8) absorbing sets with the topology shown in Fig. [6] 
for p large enough. 

Proof: The bcs for Fig. [6] has dimension 6. Using 
linearly independent cycles V\— v^— 1*3— t>i, V\— V^— V4— Vi, 
U1-U3-U4-U1, vi-va-v 6 -vi, vi-v 4 -v 5 -vi, 

V1—V5—VQ—V1, we construct the following CCM: 




Fig. 6. Configuration 4 of a (6, 8) absorbing set. This is the second of the 
two configurations with exactly one variable node that has 5 satisfied check 
nodes. 



M 



\ - 


i-2 


12 




h 







ii — 











ie - 


u 







k 






«3 - 


u 







in 























Ho - 



























is — Ho 
ig — is 







it — in 


h — «9 



(20) 



Let the top-left 3x3 submatrix of M be A. Note that A 
is exactly the CCM in ^ for the (4, 8) absorbing set of Fig. 
[T] Corollary [T] requires that dct A = mod p for the (6, 8) 
absorbing set of Fig. [6] to exist. 

As shown in Corollary [3] EAB codes always contain (4, 8) 
absorbing sets. Thus there is an assignment of EAB row group 
labels to ii, 12,13, 14, such that dct A = mod p. 

When det A = mod p, the second row of M is a linear 
combination of the first and third rows. Therefore the rank 
of M only depends on the submatrix M that contains all the 
rows of M except the second row. 

It suffices to consider the case when the labels 
i\, *2, is, U, ib, *6 adopt one of the following two assignments: 
(ii,i 2 ,iz,n,h,i&) = (x,t,w,y,z,z) or (z, t, z, y, x, w). For 
the first assignment, 



detM =(x - t)[(z - w)(i 10 - y)(i 7 - iu)(ig - is) 

- (in - z)(w - y)(i 8 - iio){h - is)], 
and for the second assignment, 

detM =(z - t)[(x - z)(iio - y){i7 - iii)(i9 - is) 

— (in — x)(z — y)(is — «io)(i7 — ig)]- 



(21) 



(22) 



For both assignments and p > 41, det M ^ mod p for 
i\ to in taking values in the set {0,1,2,3,4} such that 
the Bit Consistency constraints are satisfied. Thus, this (6, 8) 
absorbing set is not present in EAB codes for p > 41. ■ 

5) (6,8) configurations 5 and 6 - Figs, ^and^ Finally, 
we consider the case when no variable node in the absorbing 
set has all five satisfied checks. This implies one of the two 
configurations shown in Figs. [7] and [8] 

Fig. [7] contains a (4,8) absorbing set. Thus any SCM code 
avoiding the (4, 8) absorbing set will also avoid this (6,8) 
absorbing set. The EAB code, which does not avoid the (4, 8) 
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Fig. 7. Configurations 5 of a (6, 8) absorbing set. This is the first of two 
configurations with no variable nodes that have 5 satisfied check nodes. 




Fig. 8. Candidate 6 configuration of a (6, 8) absorbing set. This is the second 
of two configurations with no variable nodes that have 5 satisfied check nodes. 

absorbing set of FigJT] has 0(p 3 ) (6, 8) absorbing sets of the 
type shown in Fig. [7] for p large enough to guarantee that the 
needed (4, 8) absorbing set is present. 

Using arguments similar to those in Lemma [3] we can 
establish that for p large enough there are no (6, 8) absorbing 
sets of the type shown in Fig. [8] in either the EAB codes or 
in the SR-SCB codes. 

We omit detailed analysis of these cases for brevity. See 
ll33l for the details, which are similar to the analysis above. 

D. Summary for (4,8), (5,9) or (6,8) absorbing sets 

We have now considered all possible (4, 8), (5, 9), and (6, 8) 
configurations. The following theorem is a consequence of the 
results established above: 



TABLE I 

ROW SELECTION FUNCTIONS FOR SR-SCB CODES WITH r = 5 AND THE 
SMALLEST p PRECLUDING ALL (4,8), (5,9), AND (6,8) ABSORBING SETS. 



V 


Rate k = 1 TP=p± 


RSF 


67 


0.9263 


[0,1,2,4,17] 


73 


0.9323 


[0,1,2,3,11] 


79 


0.9373 


[0,1,2,6,7] 


83 


0.9403 


[0,1,2,3,7] 


89 


0.9443 


[0,1,2,4,11] 



Theorem 2: In r — 5 EAB codes the number of (4,8), 
(5,9) and (6,8) absorbing sets scales as 0(p 3 ). However, 
for SR-SCB codes with a properly selected RSF (such as the 
examples given in Table [E]i the relevant determinants for the 
CCMs (or submatrices of CCMs) for all (4, 8), (5, 9) or (6, 8) 
absorbing sets are nonzero allowing a p to be selected so that 
these determinants are all nonzero mod p, thus precluding 
all (4, 8), (5, 9) and (6, 8) absorbing sets. ■ 

Table [IJpresents some examples of RSFs that produce r = 5 
SR-SCB codes for which (4,8), (5,9) and (6,8) absorbing 
sets do not exist. These examples provide the smallest p that 
precludes the absorbing sets. The precluded absorbing sets can 
still appear again for larger values of p. However, in each case 
there is a p large enough that the absorbing sets are precluded 
for all larger values of p, as stated in the following corollary: 

Corollary 15: When all of the relevant determinants are 
nonzero, there is some p for which the associated absorbing 
sets are precluded for all larger values of p. ■ 
As examples, with RSF [0, 1, 2, 4, 7] and p > 179 and with 
RSF [0,1,3,4,5] and p > 271 all (4,8), (5,9) and (6,8) 
absorbing sets are precluded. 

Remark 2: For sufficiently small p's, (4, 8), (5, 9) and (6, 8) 
absorbing sets cannot be eliminated simultaneously for the 
SR-SCB codes. Here are a few examples that remove most of 
these absorbing sets for p < 50: (1) for p = 47 with RSF= 
[0, 1, 2, 3, 5], only the absorbing sets in Fig. [3] and Fig. [5]exist; 
(2) for p = 43 with RSF=[0, 1, 2, 4, 6], only the absorbing set 
in Fig. [3] exists. ■ 

E. Absorbing set spectrum in Tanner construction 

We now apply our analysis to the Tanner construction in 
||4], whose codes are in the family of SCB codes. The Tanner- 
construction codes have the parity-check matrix fZT'°i , bj with 
a, b £ GF(p) and o(a) = r and o(b) — c, where o(a) indicates 
the multiplicative order of a in GF(p). We refer to this class 
of codes as Tanner-construction codes throughout the paper. 

Lemma 5: In the Tanner graph corresponding to quasi- 
cyclic LDPC Hp°j b) bJ in [4 1, no (4,8) or (6,8) absorbing sets 
exist with parameters selected in Table I of 0. 

Proof: For the codes in Table I of (4) with girth greater 
than 6, neither (4, 8) nor (6, 8) absorbing sets exist since 
they contain length-6 cycles. For the codes in Table I of [4] 
with girth = 6, we have analyzed each code using the CCM 
approach and confirmed that in each case neither (4, 8) nor 
(6, 8) absorbing sets exist. We note that the girth = 6 codes of 
Table I in H span rates from 0.17 to 0.55. The full analysis 
is too long to include, but we provide the following example. 



14 



Consider the girth = 6 code with p=31,c = 6,r = 5as 
an example, with a = 2, b = 6. The resulting matrix H^f ^ ^ 
is a sub-matrix of the SR-SCB code with RSF [1, 2, 4, 8, 16]. 
Thus f(i,j) = d l W is achieved because the RSF enforces a 1 
and the submatrix of the SR-SCB code selects the columns 
to enforce b 3 . We set up a system of equations as before, 
and conclude that for (4, 8) absorbing sets, the only possible 
labeling for p = 31 is (ii, «2j ii, i&) = [x, t, w, y, z, z). 
Five non-isomorphic solutions, (x,y,z,w,t) = 
(2, 4, 1, 16, 8), (1, 8, 2, 4, 16), (1, 2, 4, 16, 8), (1, 4, 8, 2, 16), 
and (1, 2, 16, 8, 4) set det M = mod 31 for det M in ([12). 

Each of the solutions corresponds to the CCM equation: 



Mu 



x — t 
x — z 




t-z 
z-y 

z — w w — y 









1*3 







■■ mod 31. 



(23) 



Since detM = mod 31 and rank(M) = 2, equation (g3 
is equivalent to 



x — t t — z 
x — z 





z-y 



= mod 31. 



(24) 



If we expand the it's with ui 
equivalent to 



ji, equation d24) is 
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mod 31. 



(25) 



We denote the matrix in equation ([25) as Q n where the 
index 1 < n < 5 indexes the five realizations of Q corre- 
sponding to the five possible choices of (x, y, z, w, t) values. 
The nullspace of Q n is Q^. For any (4,8) absorbing set, 
(ill J2> jsiji) should be in U Q„. Because of the Tanner 

l<n<5 

construction, the column groups for this code are the powers 
of 6 mod 31. Denote this set by Y — {1,6,5,30,25,26}. 
The 4-tuple jz, ji) must be a selection of four distant 

elements of Y. However, since there is no vector in U Q,^ 

l<n<5 

that consists of four distinct elements of Y, (4, 8) absorbing 
sets do not exist and consequently the (6, 8) absorbing sets in 
Figs. [4] [6] and [7] also do not exist for this Tanner construction. 

Similarly, for the CCMs constructed in ( flT) for configura- 
tion 3 in Fig. [5]and for configuration 6 in Fig. [8]the null spaces 
do not include vectors that contain four distinct elements of 
Y, which precludes the absorbing sets in Fig. [5] and [8] ■ 

Codes listed in Table I of [4] are mostly moderate-rate codes 
(the rate being around 0.5). The construction of codes with 
higher rates from this code family requires more columns 
to be selected from the SCB mother matrix. However, when 
more columns are selected (by choosing a higher-order 6), the 
Tanner construction cannot guarantee the absence of certain 
absorbing sets. 

Lemma 6: In the Tanner graph corresponding to quasi- 
cyclic LDPC codes with H^'"^ in [4|, selecting a higher- 
order b to achieve higher rates can introduce (4, 8), (5, 9) and 
(6, 8) absorbing sets. 

Proof: Continuing with the example above with p — 31 and 
a = 2 which implies r = 5, we increase the rate by replacing 
6 = 6 which implies c — 6 with b = 15 which implies c = 10. 



This new choice of b yields the column selection set Y = 
{1, 15, 8, 27, 2, 30, 16, 23, 4, 29}. With this Y, (4, 8) absorbing 
sets exist because there are vectors in U that consist 

l<n<5 

of four distinct elements of Y. Similarly, (6,8) absorbing sets 
exist because corresponding solutions for Figs. [5] and [8] can be 
constructed from distinct elements of Y , ■ 

With the same row selection, the parity check matrix of 
the SR-SCB code has the parity check matrix of the Tanner 
construction as a submatrix. According to the conclusions of 
Theorem [2j we can modify Tanner construction by using an 
optimized row selection that precludes all the (4, 8), (5, 9) and 
(6, 8) absorbing sets for appropriate values of p. 

Lemma 7: In the Tanner graph corresponding to quasi- 
cyclic LDPC codes (in [4|) described by H^'°l b y with higher 
rates and large enough p, all of (4, 8), (5, 9) and (6, 8) ab- 
sorbing sets can be eliminated by only modifying the mapping 
sequence of the rows. 

Proof: By replacing the row mapping sequence g r (i) = 
a" 1 - 1 by g r (i) = [0,1, 2, 4, 17] and choosing p = 67, the quasi- 
cyclic LDPC H 5 '°S h) s becomes a new class of SCB LDPC 

codes with f(i, j) — q r (i) ■ V . Since H 5 ' .^ is a submatrix 
of the mother matrix H^'?., it follows that (4, 8), (5, 9) and 
(6, 8) absorbing sets can be structurally eliminated for p large 
enough using the exactly the same analysis as in Sections [IV-A| 



through IV-C 



Remark 3: We already showed that one can easily check the 
existence of certain absorbing sets in the Tanner construction 
101 (Lemmas [6] and [7). Other popular quasi-cyclic LDPC 
codes, such as codes in [5| and J6) can also be viewed as 
being constructed as a selection of certain rows and columns 
of the SCB matrix. Codes in Section III-B2/B3 of use 
the SCB structure with parity-check matrix H r ' c f ,. where 
f(hj) = £ " h < i < r — 1, < j < c — 1 or 
f (ij) = (a^-lXV-l), < i < r-l, < j < c-1. Array 
codes [6 1 use the first r row groups of H^'f.j- The approach 
developed for the SCB construction can therefore be easily 
applied to codes in J5] and J6). 

F. Equivalence Classes for SR codes 

We now present three code equivalence conditions and use 
them to fully characterize the p = 7, r = 5 case. 

Since the order of the elements in the RSF only permutes 
the rows of the matrix H r ' c r , and does not change the code 
properties, we can assume that the RSF vector [ai, a r ] 
is ordered in ascending order. 

Consider a difference matrix D of the RSF, where D i3 = 
a,j — ^ mod p, 1 < i, j < r. If D = D or D is the reflection 
of D on its antidiagonal, we say that D and D are equivalent 
difference matrices. The following lemma establishes some 
absorbing-set equivalence classes. 

Lemma 8: The following pairs of RSFs are equivalent in 
the sense that they have the same absorbing sets: 

1) [a%, a,2, a r ] = [oi, 02, a r ] + constant mod p 

2) [a%, 02, a r ] = [a±, 02, a r ] x constant^O mod p 

3) [a%, a,2, a r ] = [5i, 02, a r ] if they have equivalent 
difference matrices. 
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TABLE II 

The 63 class leaders (excluding and 1) for RSF equivalence 

CLASSES PRECLUDING (4,8), (5,9) & (6,8) ABSORBING SETS FOR p = 67. 



2,4,17 


2,12,17 


2,16,48 


3,8,42 


3,16,39 


3,27,49 


3,54,63 


2,4,48 


2,12,24 


2,17,45 


3,8,54 


3,17,30 


3,27,52 


4,13,47 


2,4,53 


2,12,38 


2,24,27 


3,8,57 


3,17,47 


3,29,50 


4,13,61 


2,4,64 


2,12,41 


2,25,41 


3,10,50 


3,17,49 


3,29,51 


4,16,19 


2,5,17 


2,13,42 


2,28,38 


3,10,54 


3,18,39 


3,30,38 


4,19,42 


2,5,44 


2,13,44 


2,31,32 


3,15,18 


3,20,21 


3,30,40 


4,20,26 


2,5,60 


2,13,48 


3,4,17 


3,15,46 


3,21,40 


3,42,49 


4,44,54 


2,9,32 


2,16,32 


3,4,20 


3,15,63 


3,21,59 


3,46,50 


4,47,48 


2,9,44 


2,16,38 


3,8,37 


3,16,20 


3,27,29 


3,52,59 


4,53,61 



Proof: the existence of any particular absorbing set hinges 
on whether the determinant of the associated CCM (or certain 
submatrices of the CCM) is zero. Since the determinant is 
only a function of the differences between the elements of the 
RSF, if two RSFs [ai, a®, • a r ] an d [Si, 5,2, ■ 5 r ] share any 
of the three equivalent conditions, [aj., <i2> •••) <v] leads to a 
zero determinant if and only if [ai, a^, a r ] leads to a zero 
determinant. ■ 

Remark 4: Since the null space of CCM also only depends 
on the difference of the column-group indices, analogous 
equivalence conditions can be established column-wise. ■ 

The following is a consequence of Lemma [8] 

Corollary 16: Any RSF is equivalent to an RSF 
[0, 1, x, y, z], where x,y,z are positive integers. 

Proof: With condition ( 1 ) in Lemma [8] a mapping vector 
is equivalent to any mapping vector that differs by a constant 
in GF(p). Thus we can subtract the smallest value in the 
mapping vector and obtain a in the mapping vector. Then, 
with condition (2) in Lemma [8] and the multiplicative property 
of GF(p), we can multiply the mapping vector by some value 
in GF(p) to make one of the non-zero elements equal to 
1. Therefore any mapping vector is equivalent to a mapping 
vector that contains and 1. ■ 

This result enables a reduced search of structured matri- 
ces with good error-floor properties. For example, a row- 
mapping vector [0, 1, 2, 4, 6] is equivalent to [1, 2, 3, 5, 7], and 
to [0,2,4,8, 12]. 

Here is an example for p = 67. There are ( 6 5 7 ) = 9657648 
RSFs, but there are only 2192 equivalence classes. These 
equivalence classes can be represented by their class leaders 



which start with and 1, according to Corollar) 16 Only 63 
out of 2192 classes preclude (4,8), (5,9) and (6,8) absorbing 
sets. Table [ill lists all the 63 class leaders, where the common 
and 1 are omitted to save space. To generate all the equivalent 
RSFs, we can couple the and 1 to each 3-element set and 
form a 5-element RSF vector. Then by multiplying i and/or 
adding j in GF(p), where i,j — 1,2, ...,p — 1, we can 
reconstitute all the equivalent RSF vectors that preclude (4,8), 
(5,9) and (6,8) absorbing sets. 

V. Theoretical results for r = 4 

This section provides an example with r = 4 (four row 
groups) that shows how to design an SCB code with a specified 
circulant matrix that eliminates the dominant absorbing sets by 
selecting rows and columns from the SCB mother matrix to 



force the CCMs associated with the dominant absorbing sets 
to have zero-dimension nullspaces. 

In addition to the EAB codes and SR-SCB codes introduced 



in Section IV removing a few column groups from an SR- 
SCB code provides further improvement. Hence, shortened SR 
(SSR) SCB codes form our third class of SCB codes. The 
parity-check matrix for these codes is H^'j, with f(i,j) = 
g r (i)'9c(j) where g c (j) is called the column-selection function 
(CSF). Note that for a pxp circulant matrix, EAB and SR-SCB 
codes have p column groups (p 2 binary columns), but SSR- 
SCB codes have fewer column groups since g c (j) selects a 
subset of the possible column groups. 



Section V-A identifies the (6,4) absorbing sets as dominant 
for EAB codes with r = 4. Section |V-B| analyzes the three 
possible (6, 4) absorbing set configurations and shows how 
carefully selecting four row groups from the SCB mother 
matrix can eliminate two of the three possible configurations. 
Section [V-C| provides an efficient provable algorithm to elim- 
inate all (6, 4) absorbing sets by combining the row selection 
of Section IV-BI with column selection in which some column 
groups of the SCB mother matrix are removed. 

A. Identifying the dominant absorbing sets 

From the previous results in |[T5l . (6,4) absorbing sets are 
the smallest possible structure for EAB codes with r — 4 for 
p > 19. Hardware simulations |3j also demonstrate that (6, 4) 
absorbing sets are the dominant cause of the error floor for 
example r = 4 EAB codes. 

Based on these results, a key goal will be to design an r = 4 
SCB code that avoids all (6,4) absorbing sets. The lemma 
below establishes that the new code design approach does not 
introduce other smaller absorbing sets that were avoided by 
the codes discussed above. 

Lemma 9: In the Tanner graph corresponding to 
there is no absorbing set smaller than (6, 4) for p large enough 
with a careful choice of the row-selection function g r (i). 

Proof: The smallest possible absorbing sets for an Hp'g r^.j 
SCB code are (4,4), (5,2), (5,4) and (6,2) absorbing" sets 
(cf. El). The (4, 4), (5, 4) and (6, 2) absorbing sets for r = 4 
are sub-graphs of the (4, 8), (5, 9) and (6, 8) absorbing sets, 
respectively, for r = 5. From the analysis of absorbing sets for 
r = 5 in the previous section (see also ifTOl ). a careful choice 
of the r = 5 row-selection function (RSF) a(i) eliminates 
the (4, 8), (5, 9) and (6, 8) absorbing sets for p large enough. 
Taking any 4-element subset g r (i) of such an 5(i), for example 
the RSF g r (i) where \g r (0),g r (l),g r (2),g r (3)] = [0,1,3,4], 
as a subset of [0, 1, 3, 4, 5] that can eliminate the (4, 8), (5, 9) 
and (6,8) absorbing sets for p > 271, avoids the (4,4), (5,4) 
and (6, 2) absorbing sets for r = 4. Any RSF that avoids the 
(4,4) absorbing set also avoids the (5,2) absorbing set |15|. 
Thus the resulting E^,.,. avoids (4,4), (5,2), (5,4) and 
(6, 2) absorbing sets. ■ 

Remark 5: Since the SSR code is a shortened version of 
the SR code, obtained by removing certain variable nodes, no 
smaller absorbing set will be introduced in the SSR code. ■ 



In Section V-B we show that SR codes always have (6, 4) 



absorbing sets, irrespective of the choice of g r (i). Avoidance 
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Fig. 9. Depiction of the first (6, 4) absorbing set configuration. 

of all such configurations using shortening is the subject of 
Section IV-Cl 



B. (6, 4) absorbing sets in SR-SCB codes 

Three distinct configurations of (6, 4) absorbing sets are 
possible for r — 4 SCB codes. This section shows which of 
these configurations are possible in EAB and SR-SCB codes. 
The first configuration exists in the EAB code and in every 
possible SR-SCB code. The second configuration exists in the 
EAB code but can be avoided by a proper choice of the RSF 
for the SR-SCB code. The third configuration does not exist 
in either the EAB code or the SR codes. 

1) The first (6,4) configuration: Fig. [9] shows the first 
configuration of a (6, 4) absorbing set in an r = 4 SCB code. 
The following lemma establishes that the EAB code and all 
SR codes have this configuration. 

Lemma 10: In the Tanner graph corresponding to the EAB 

p,/(v) 



code and all SR codes with H*'^,, ^ there are (6, 4) absorbing 



sets for any p with the configuration shown in Fig. [9] 

Proof: Using the technique of Section III-B| we con- 
struct the CCM for this configuration. The bcs for 
Fig. [9] has dimension 5. Using the five linearly indepen- 
dent cycles: v\— v 2 — v^— V\, V\— v 2 — ^4— «i, V\— v 2 — v$— v\, 
V\— V3— Vq— V4— V\, vi—v^—vq—v^—vi, we construct the fol- 
lowing CCM: 



M 



14 — is ib — ii 

24 — ie ie — ii 

24 — 17 17 — ii 

ii — is is — «2 is — ig 

ii — ii, iio — iz *9 — iio 



(26) 



To avoid this absorbing set we need to force the determinant 
of the CCM to be nonzero. Note that det(M) is computed as 

M11M23M34M42M55 

- Mi 2 (M23M3iM 45 M54 - M 2 iM3 4 (M 43 M55 - M 45 M 5 3)), 

where My denotes the entry in M. 

From the proof of Lemma 8 in ifTSIl . there are 
only two non-isomorphic row-group labelings for the 
check nodes of Fig. [9] It is sufficient to consider only 

{h-,i2,h,h-,h,k,i7,i&,k,ixo) = (x,y,z,w,y,z,x,z,x,y) 
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Fig. 10. Depiction of the second candidate (6, 4) absorbing set. 

or (x, y, z, w, y, z, x, z, w, y) with x, y, z, and w taking distinct 
values in the range of RSF. The first labeling yields 



M 



w 

w ■ 

w - 



V 

z 

X 



y-x 



x — z 






z-y 


x - y 

y-x 






x — z 


y - z 







z — x 
x-y 



(27) 



~u 2 ~ 




113 












_ue_ 





and det(M) = regardless of the specific values of w, x, y, z 
taken as mutually distinct integers mod p. 

Thus there exists a non-zero solution to M • u = mod p. 
One solution to this equation is 

(x - y)(z - y)(x - z) 
(w - y)(z - y)(x - z) 
(w - z)(y - x)(x - z) 
(w - x)(y - x)(z - y) 
(w - y){x - z)(z -y) + (y- x)(w - z)(x - y)_ 

(28) 

For this absorbing set, the Check Consistency condition 
requires U2 ^ 0, M3 ^ 0, U4 ^ 0, U5 ^ 0, ^ 113, 
U2 7^ U4, U2 7^ U5, U3 7^ Uq, U4 7^ uq, and U5 ^ uq. These 
requirements as well as the Bit Consistency inequalities are 
met since x, y, z, w are mutually distinct. 

The solution in ( f28] i satisfies the Bit, Check and Cycle 
Consistency constraints. This is a sufficient condition for the 
existence of a (6, 4) absorbing sets with the configuration of 
Fig. [9] Any four distinct values between and p — 1 for 
{x, y, z, w} identify a labeling of this first type that identifies 
an absorbing set in the EAB code and every SR SCB code. 

Consider . . . ,i w ) = (x,y,z,w,y,z,x,z,w,y), the 

second labeling. In this case det(M) 7^ for {x,y, z,w} = 
{0, 1, 2, 3}. Thus det(M) =k mod p for p large enough, and 
there is no such (6, 4) configuration with this labeling in the 
EAB code for p large enough. The EAB code is one possible 
SR code. Other careful choices of the RSF, produce other SR 
codes that likewise do not have this (6, 4) configuration. ■ 

2) Two additional configurations: Fig. [T0| shows the second 
possible configuration of a (6,4) absorbing set in an r = 4 
SCB code. Similar arguments to those above establish that the 
EAB code has this configuration but well-designed SR codes 
avoid it. One such example is the RSF [0, 1, 3, 4] which avoids 
this configuration for p > 31. 

Fig. 11 shows the third configuration of a (6,4) absorbing 
set in an r = 4 SCB code. Similar arguments to those above 
establish that neither the EAB nor well-designed SR codes 
have this configuration for p sufficiently large. In fact, p > 13 
is sufficient for the absence of this configuration in the EAB 
codes. SR-SCB codes can avoid the configuration in Fig. 11 
with a careful choice of the row mapping. One such example is 
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Fig. 11. Depiction of the third candidate (6, 4) absorbing set. 

the EAB RSF itself. Another is to choose the RSF [0, 1, 3, 4], 
which can avoid this configuration for p > 19. 



Algorithm 1 Greedy column-cutting algorithm. 
1: % Initialization: C begins with all columns. 

2: C = {0,l,...,p-1}. 

3: Collect all the vectors \j\,j2---ije\> that form the config- 
uration in Fig. [9] with j n G C, n = 1, 6 and they form 
a set W= V. 

4: % Proceed one column at a time removing columns that 
% preclude vectors in V until no vector in V is possible. 

5: while \W\ > do 

6: Find the most frequent j in W, say j m . 
7: Replace C by C \ j m . 

8: Remove every \jx,j2»-, je] that involves j m from W. 
9-. end while 

10: % C contains the column groups of the designed code. 



C. Eliminating (6,4) absorbing sets with shortening 
For a sufficiently large p, well-designed SR codes avoid the 



Algorithm 2 Column-adding algorithm. 
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How- 



(6, 4) absorbing set configurations in Figs. [TO] and 
ever, as shown in Lemma [10] SR codes cannot eliminate the 
(6, 4) absorbing set configuration in Fig. [9] We now consider 
shortened SR (SSR) codes that retain only certain column 
groups from the SCB mother matrix (reducing the rate). A 
well-chosen column selection g c (j) allows the Tanner graph 



corresponding to H 



sets. 



P,9r(i)-ffc(i) 



to avoid all (6, 4) absorbing 



We begin with an SR code using well-selected RSF g r (i), 
for instance [0, 1, 3, 4], that already avoids the (6, 4) absorbing 



set configurations in Figs. [10] and 1 1 for p large enough. We 
then choose a CSF g c (j) to also avoid the (6, 4) absorbing set 
configurations in Fig. [9] Choosing a column selection g c {j) 



reduces to choosing a submatrix of H 



4,p 

P,9r-(i)j 



by eliminating 



certain variable nodes. This operation cannot introduce smaller 
absorbing sets. 

One solution to M-u = mod p is equation p8) . The rank 
of M in |27) is 4, and therefore this single solution forms a 
basis of the null space. Multiplying u by a constant c, for 1 < 
c < p— 1, also results in a solution. These p— 1 solutions cover 
all of the feasible solutions described by the null space. For 
fixed u\, M5, we can choose j\ from 0, 1, ...,p—l and obtain 
J2, j6- Thus there are p(p — 1) ways to find j\ to j$ for a 
fixed {x, y, z, w}. Since there are 4! ways to assign check node 
labels based on the set {x, y, z, w} for a fixed row mapping, 
there are at most 24p(p— 1) possible vectors \ji,j2--ije] that 
can form the configuration in Fig. [9] These vectors form the 
set V of vectors, which completely characterizes the instances 
of this absorbing set configuration. 

If a CSF is applied, each variable node group label j is in 
a set J where J C {0, 1, ...,p — 1} and we can only choose 
Uuh—Je] such that j m e J,m = 1,...,6. There are ('g 1 ) 
possible [ji, j2--., je] vectors and they form a set of vectors 
V. If V PI V = 0, the new code does not have the (6, 4) 
configuration in Fig. [9] We can find such a CSF with the 
greedy column-cutting procedure described in Algorithm [T] or 
the column-adding procedure described in Algorithm [2] 

Remark 6: A similar technique could be applied to increase 
the girth [28 1 instead of eliminating the smallest absorbing set. 



1: % Initialization: Select an initial six columns for C that 

% do not form a vector in V. 
2: J = {0,l,...,p-1}. 

3: Collect all the ordered vectors j2---, je], that form the 
configuration in Fig. [9] with j n 6 J,n = 1, 6, into the 
set V. 

4: Choose a distinct 6-element set C, C C J. 

5: while for some ordering [ji, j2---, je] of the elements of 

C, fiuh-M € V do 
6: Choose another distinct 6-element set C randomly such 

that C C J. 
7: end while 

8: % Proceed one column group at a time, adding columns 
% to C that do not introduce the absorbing set. 

9: J = {0,l,...,p-l}\C 
10: while |J| > do 
11: Select a j m randomly from J. 

12: if every $1,32-, h] & V for every j 1: j 2 ..., Je G {C U 

j m } then 
13: C = C U j m 
14: end if 

15: J — J\ j m 

16: end while 

17: % C contains the column groups of the designed code. 



However, simply increasing the girth would not guarantee a 
better performance, see e.g., [29 1. ■ 



D. SSR codes with Tanner construction 

This section shows that for r = 4, the row-selection function 
g r (i) of the Tanner construction J4] will always introduce 
(4, 4) absorbing sets for the case set forth in the following 
lemma. For a an element of GF(p) let o(a) denote its 
multiplicative order in GF(p). 

Lemma 11: In the Tanner graph corresponding to the quasi- 
cyclic LDPC code with the parity-check matrix Hp'g~(iy g uy 
where g r (i) = a l ,o(a) = 4,g c (j) = V ,o(b) = p - 1, (4,4) 
absorbing sets as shown in Fig. [12] always exist. 
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(j..k.) 



Fig. 12. Depiction of the (4,4) absorbing set configuration. 

Proof: The proof uses the same techniques as that of 
Corollary [2] Details are omitted for brevity. ■ 

Remark 7: Tanner constructions with o(b) < p — 1 and 
r = 4 sometimes can avoid all (4, 4) absorbing sets by 
either increasing the girth or by making the intersection 
between null space of M and variable node space to be 
empty. However, (4, 4) absorbing sets may still exist in some 
constructions with o(b) < p — 1 and r = 4. For example, with 
p = m,g r {i) = H\g c {]) = 5*, o(ll) = 4,o(5) = 30, the 
resulting code has (4, 4) absorbing sets. 

VI. Discussion for r = 3 

We briefly remark on the case where the column weight is 
3. First, we note that for the EAB codes with r — 3, the (3, 3) 



and (4. 2) absorbing sets in Fig. 13 are the smallest ones, 031 



It is easy to show that the same absorbing sets are also the 
smallest possible for SR-SCB codes (of girth 6). The CCM 
expression describing a (3, 3) absorbing set now takes on a 
particularly simple form: 

h(j3 ~ ji) + h{j2 - h) + - h) = mod p, 

which always has the solution j\ = 13, — h an d 
J3 = h- The number of (3,3) absorbing sets scales as 0(p 3 ). 
This expression precisely describes a length-6 cycle and thus 
establishes a one-to-one relationship between (3, 3) absorbing 
sets and length-6 cycles. It is therefore necessary to increase 
the girth for the elimination of (3, 3) absorbing sets. The 
resulting shortened codes would have much lower rates and 
as such are not of interest in this paper. 

Similarly, we can apply the CCM approach to the (4, 2) 
absorbing sets and easily prove their existence in any EAB 
and SR-SCB codes. The CCM equation now becomes 



M 



h -h 12- k 
15 - h h — u 



(29) 



and always has a nontrivial nullspace. Thus this absorbing set 
cannot be avoided in the SR-SCB codes. The number of (4, 2) 
absorbing sets scales as Q(p 3 ). 

VII. Simulation Results 

In this section we experimentally demonstrate the perfor- 
mance improvement achieved by the proposed modifications. 
We show that it is consistently valid across various choices of 
decoding implementations and platforms both in software and 
in hardware. Simulations are performed over additive white 
Gaussian noise (AWGN) channels. The limited-precision de- 
coders in the simulations use 200 iterations and a Q4.2 
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Fig. 13. Depiction of the (3, 3) and (4, 2) absorbing set configurations. 
TABLE III 

Software error profiles for three (2209, 1978) codes: the EAB 

CODE, AN SR-SCB CODE AND HUANG'S CONSTRUCTION. THE TOTAL 
NUMBER OF COLLECTED ERRORS IS DENOTED N.E. THE NUMBER OF RUNS 
IS DENOTED N.R. 







V = 


47 (2209, 1978) EAB code 






SNR 


n.r. 


n.e. 


(4,8) 


(5,9) 


(6,8) 


(6,10) 


(7,9) 


(7,11) 


(8,6) 


5.2dB 


3.3E8 


378 


17 





174 





68 








5.4dB 


3.8E8 


230 


5 





113 





72 












V = 


47 (2209, 1978) Huang code 






SNR 


n.r. 


n.e. 


(4,8) 


(5,9) 


(6,8) 


(6,10) 


(7,9) 


(7,11) 


(8,6) 


5.2dB 


3.0E8 


350 


51 


25 


106 


6 


70 


11 


3 


5.4dB 


5.2E8 


301 


39 


29 


115 


2 


61 


15 







p = 47 (2209, 1978) SR-SCB code RSF = 


[0,1,3,8,19] 




SNR 


n.r. 


n.e. 


(4,8) 


(5,9) 


(6,8) 


(6,10) 


(7,9) 


(7,11) 


(8,6) 


5.2dB 


8.1E8 


241 








4 


69 


9 


13 


29 


5.4dB 


2.1E9 


165 











56 


2 


15 


8 



quantization, 4 bits to the left of the radix point to represent 
integer values and 2 bits to the right of the radix point to 
represent fractional values. We simulate the sum-product al- 
gorithm (SPA) and the soft xor algorithm (SXOR) [34] for 
different codes. We also identify the absorbing sets responsible 
for each error in the error floor and report an error profile for 
each simulation. As expected, when an SR-SCB code or SSR- 
SCB code is designed to preclude certain absorbing sets, they 
do not appear in the error profile obtained by simulation. 

A. SCB codes for r = 5 

Consistent with the results presented at the end of Section 
IV^Al Table Inl] and Fig. 



14 show the software, limited- 



precision SPA performance of two (2209, 1978) SCB codes: 
the EAB code and an SR-SCB code. Both codes have the 
following parameters: check node degree = 47, and bit node 



(2209,1 978) LDPC AWGN fixed point Q4.2 




Fig. 14. Comparison of the (2209, 1978), r = 5, rate-0.895 EAB, SR-SCB 
LDPC codes and Huang's construction [25 ] using SPA decoder. 
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degree = 5. The EAB code uses g r (i) = i, and SR-SCB code 
uses the RSF [0, 1, 3, 8, 19]. Another (2209, 1978) quasi-cyclic 
code using Huang's construction 11251 with the same code 
parameters as our SCB codes is constructed and simulated 
with the same decoder. We also simulated all three codes using 
a limited-precision SXOR decoder in software. SXOR results 
are omitted in the interest of space. 

Irrespective of the decoding algorithm (SPA or SXOR), the 
SR-SCB code that provably eliminates all (4, 8) absorbing 
sets performs the best among the three codes at every SNR 
point measured. Both EAB and Huang's code contain (4, 8) 
absorbing sets. At measured SNR points, Huang's code is 
always in between the EAB and SR-SCB codes in terms 
of the performance. Huang's code offers visible performance 
advantage over the EAB code under SXOR, in part due to the 
suppression of the effect that the (4, 8) absorbing sets of this 
code have under SXOR decoding. Under SPA (and as shown 



in Figure 14 1 this performance advantage vanishes. 

In addition to software simulations, we also performed SPA 
FPGA simulations for (2209, 1978) EAB and SR-SCB codes 
(results omitted in the interest of space, see also [10]). These 
results also demonstrate that the BER improvement is about 
one order of magnitude at 5.8 dB for SR-SCB codes relative to 
EAB codes. This is the SNR point where SR-SCB codes reach 
BER of 10~ n . The majority of collected errors for the EAB 
code in the low BER region is again due to (4, 8) absorbing 
sets. The FPGA simulations collected more errors overall so 
that (5, 9) absorbing set errors were observed with both the 
EAB and SR-SCB codes. 



While Table III and Fig. 14 showed that performance 



improvement can be obtained by precluding only (4,8) ab- 



sorbing sets, Table IV and Fig. 15 show hardware simulations 
demonstrating the benefit of precluding all (4, 8), (5, 9) and 
(6, 8) absorbing sets. Here we compare longer block-length 
EAB and SR-SCB codes wither r = 5 and p = 67. The SR- 
SCB code precludes all (4,8), (5,9) and (6,8) absorbing sets 
using the RSF [0, 1, 2, 4, 17] with p = 67 which was identified 
in Table [I] of Section IV-D| as completely avoiding these three 
absorbing sets. As expected, the error profile as shown in 
Table |IV] shows that the SR-SCB code completely eliminates 
the (4,8), (5,9) and (6,8) absorbing sets. Fig. 15 shows the 
corresponding BER improvement. 

B. SCB codes for r — 4 

Fig. [16] and Table [V] show the performance of (2209, 2024), 



p = 47 EAB code and SR-SCB codes both with check 
node degree = 47 and bit node degree = 4. The EAB code 
has g r (i) = i and the SR-SCB code uses RSF [0,1,2,4]. 
Consistent with the theoretical analysis, the (6, 4) absorbing 
sets dominate the error floor of the EAB code. Precluding 



the (6,4) absorbing sets of Figs. 10 and 11 provides the SR 
performance improvement. Precluding the (6,4) absorbing set 
of Fig. [9] as well requires shortening which lowers the code 
rate somewhat. Such a code is discussed in the next paragraph. 
16 also compares the performance of a high-rate 



Fig. 



TABLE IV 

Error profiles for the EAB SCB (4489, 4158), code (top), and 

THE SR-SCB CODE (BOTTOM), BOTH WITH p = 67. N.E. IS THE NUMBER 
OF COLLECTED ERRORS. N.R. IS THE NUMBER OF RUNS. 



(4489, 4158) EAB code 



SNR 


n.r. 


n.e. 


(4,8) 


(5,9) 


(6,8) 


(6,10) 


(7,9) 


(7,11) 


(8,6) 


(8,8) 


5.6dB 


1.1E9 


150 


67 


17 


22 


7 


6 


5 


6 


6 


5.8dB 


2.1E9 


139 


83 


18 


16 


6 


5 


1 


3 


3 


6.0dB 


4.3E9 


131 


77 


18 


22 


5 


1 


1 


2 


1 


6.2dB 


8.6E9 


107 


85 


10 


5 


4 


2 











(4489, 4158) SR-SCB code RSF = 
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LDPC (p=67,r=5) AWGN fixed point Q4.2 



-N=4489 EAB code ratc-0.9263 
-N=4489 SR code rate=0.9263 



quasi-cyclic (QC) code under the Tanner construction QJ 
with a similar-rate shortened SR (SSR) code that precludes 



Fig. 15. Performance comparison of (4489, 4158) EAB and SR-SCB codes 
using limited-precision SPA decoder. 



all (6,4) absorbing sets. The QC code has the following 
parameters: p = 61, f(i,j) = a 1 -V , a — 11, 6 = 5, o(a) = 4 
and o(b) = 30, where o(a) indicates the multiplicative order 
of a in GF(p). Using the CCM based analysis, one can show 
that this code does not have (6, 4) absorbing sets, although 
it does have (4, 4) absorbing sets (due to an inappropriate 
row mapping). The SSR-SCB code is obtained by shortening 
the SR-SCB code with parameters p = 79 and RSF = 
[0,1,3,4] using Algorithm 1. This SSR-SCB code has CSF 
[2, 6, 7, 14, 17, 18, 22, 26, 27, 30, 36, 37, 38, 46, 47, 49, 55, 56, 
57,58,61,62,65,66,67,76,77,78]. We thus obtain a code 
with the same variable-node degree as the QC code, a similar 
block length (N = 2212 vs. N = 1830 for the QC code) 
and a similar rate (0.8585 vs. 0.8683 for the QC code). The 
SSR-SCB code provably eliminates (6,4) absorbing sets 
without introducing smaller absorbing sets. The profiles in 
Table [V] confirm this claim. Similar analysis can be applied 
to the codes in J5) as well. 

The (1944, 1620) quasi-cyclic code from 802.1 In standard 
is also simulated as a reference. This code has larger block 
length and lower rate than our SSR-SCB codes. The error floor 
is dominated by the small absorbing sets such as (3,3) and 
(4, 1) absorbing sets and the (2, 2) trapping set which consists 
of a pair of degree-2 variable nodes. This code has a relatively 
early error floor compared to the SSR code. 

Lastly, we plot the same SSR, Tanner-construction QC, and 



802.1 In codes again in Fig. 17 now with a full precision 
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TABLE V 

Error profiles for (2209, 2024) EAB (p = 47), SR-SCB (p = 47), 
AND SSR-SCB (p = 79) CODES. 

SNR | n.r. | n.e. | (6,4)| (6,6)| (7,4)| (8,2)| (8,4)| (9,4)| (10,4) (12,4) 
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is a deterministic method that can provably eliminate certain 
absorbing sets in a large family of circulant-based codes, 
and can do so without changing the code's properties such 
as girth, rate and implementation complexity. This approach 
thus offers a class of codes with provably better performance 
than some known constructions. Theoretical findings were 
substantiated by experimental results showing consistent per- 
formance improvement over a range of decoding algorithms 
and implementation platforms. An interesting problem for 
future investigation would be to apply the CCM method to 
other graphical structures of interest, e.g., trapping sets. 



LDPC ;r 4) A WON I xo til 04 2 
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«-N 2209 SR code, rate 9163 
s-N 1830 Tanner OC code rale I 
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* -N 2212 SSR code rale 8585 
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Fig. 16. Performance comparison of EAB, SR, SSR codes, a code from |4|, 
and a code from the 802.1 In standard with a decoder using 4.2 limited- 
precision SXOR decoding. 

LDPC (r=4) AWGN full precision 



-N 1830 Tanner QC code soil, rate 8 
-N 2212 SSR code soil, rate 8585 
-N 1944 802.11nsoft.rate 8333 



Fig. 17. Performance comparison of the SSR code, the code from |4| and 
the code from the 802.1 In standard as in Fig. [16] but with a full-precision 
floating-point SPA decoder. 



decoder. We notice that (as expected) the performance of full 
precision-decoded codes improves for all codes considered. It 
is interesting to observe that the SSR code again eventually 
outperforms the 802.1 In code, although the improvement is 
much less than in the limited precision case, as the importance 
of absorbing sets is not as significant under full precision as 
it is with a limited-precision decoder. 

VIII. Conclusion 

We introduced a novel cycle consistency matrix (CCM) 
description of dominant absorbing sets to guide code design 
and analysis of circulant-based LDPC codes. Our approach 
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